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FOREWORD 


This report has been prepared for NASA Manned 
Spacecraft Center under Contract NAS9-12578. The 
technical monitor of this effort has been Claude A 
Graves of the MSC Mission Planning and Analysis 
Division. 

The principal investigator for this effort 
has been Raymond Morth. Neal Carlson, Richard Ku, 
and William Widnall consulted on the navigation 
study. 
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CHAPTER ‘ 1 


INTRODUCTION AND SUMMARY 


1.1 Study Objectives 

This report is concerned with the steering of the 
Shuttle vehicle in the final several hundred miles prior 
to landing. The steering starts just prior to the 
transition maneuver from high to low angle-of -attack and 
ends with the touch-down of the vehicle on the landing 

strip. 

The objective of this study is to extend and complete 
the analytic design that was made as part of the Phase 
effort [1] . Here simple formulae were used to represent 
the vehicle performance characteristics. These formu ae 
were used in a guidance scheme that is characterised by a 
spiral turn to dissipate the excess potential e ^ er 5Y 
(altitude) prior to a standard straight-m final app 
That study was fairly complete but several issues that 
were not resolved due to lack of time are the subgect 
of this report. 


Most important of these improvements and extensions 

are: 

3 _ % Tailor the steering with a view to pilot desires 
and capabilities 


2. Add a flare maneuver logic 


3. Design steering logic for transition phase 

4. Incorporate phugoid damping 

5. Decouple the center-of -curvature control from 
ranging control 

6. Incorporate the use of drag brakes 

As a second part of this study, this improved logic is 
to be used to define targeting procedures for the nominal 
Shuttle trajectory. 

Thirdly, the logic is to be stressed by both variation 
in key parameters and by realistic system errors. These 
include: 
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Aerodynamic parameters 

2. Atmospheric variations 

3. Winds (both steady and gusts) 

4. Initial condition variations 

•The error performance is to be characteristic of a navigation 
system which includes an inertial measurement unit and an 
on-board navigation computer. This on-board system is to 
be updated by both conventional landing aids, i.e., 

1. barometric and radar .altimeters 

2. VOR and DME of a TACAN type system 

3. ILS system, both glide slope and localizer 

and also precision DME equipment as characterized by 
the Cubic CR 100. 

1.2 Key Design Considerations 

Several design considerations were listed at the 
start of the study. They are for the most part met 
naturally by the analytic nature of the guidance design. 
These considerations as numerated by the contract monitor 
are: 

1. Must be capable of using all the vehicle 
ranging ability 

2. Compatible with manual procedures 

3. Compatible with ground control procedures and 
monitoring 

4. Capable of converging in final approach target 

5. Maintain an energy reserve for the nominal case 

6. Maintain the ability to compensate for winds 

7. Use minimum navigation aids in contingency cases. 

8. Targeting flexible to permit landing at alternate 
runways 
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Capable of performing category III landings 

10. Maintain subsonic velocities 

11. Constrained to alleviate over-pressure during 
transitions 


1.3 Background, Other Design Approaches 

\ 

There have been a large number of guidance schemes 
proposed for approach and landing phase of the Space 
Shuttle, and for all other guidance phases for that 
matter. In fact, there are so many that it will be 
difficult to select the primary system from all the 
candidates. On the other hand, this decision may not be 
too difficult because all the guidance schemes "work" and 
any one that is used will do the job. It is only in the 
way the different schemes will handle the many contraints 
and adapt to system errors that one system may show 
superior performance over the others. 

The many schemes are roughly divided into two 
categories; non-turning and turning approach patterns. 

The Bell system [ 2 ] is characterized by a single turn 
on to final approach. Detailed calculations based on 
poorly defined system parameters are characteristic of 
this system. The Draper Lab [ 3 ] defines a system with no 
turns relying on drag brakes to dissipate the excess energy. 

If the drag brakes do not have sufficient power such a mode 
is not possible. But direct control of the kinetic 
plus potential energy as proposed in this system seems to 
have some merit, particularly in the presence of large 
disturbing winds. 

The turning approach can further be divided into 
three categories. There is a two— turn or so called 
"racetrack" system as proposed by several investigators 
at MSC [4,5]. Here there is a turn on to a down -wind 
leg away from the runway until sufficient altitude is 
lost so that the second turn on to final approach is 
appropriate. Then there is the so called cone cylinder 
approach of MACDAC [6], in which the vehicle dissipates the 
energy by first flying about a large radius circle, 
then at the appropriate time changes the flight to a small 
circle tangent to the final approach path. 

Finally, there is the spiral descent path which is 
the subject of this report and also of earlier work by 
Sperry Rand [7] . The difference in the two spiral guidance 
techniques is that this approach is characterized by 
simple analytic formulae which are used to represent the 
vehicle characteristics and the guidance loop is then 
closed using these formulae. It is of interest to note 
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that the first flight tests of the Sperry Rand system look 
very similar to the two-turn system [8]. It is also 
significant to note that this system has demonstrated 
automatic landings to touchdown on a real aircraft 
{a CV 990 modifi.ed to fly like the Space Shuttle) using 
real equipment {an inertial measurement unit (IMU), an 
airborne computer, a radar altimeter, and other conventional 
landing aids) . 

In the original work for this study [1] , the North 
American .vehicle had a much lower wing loading. With the 
resultant lower velocity at altitude, sometimes up to 5 
turns would be required. Here, as with similar studies 
on the Grumman Phase B vehicle, no more than two turns 
is observed and it is usually one. But should this capability 
return {it never does as projects evolve) the present logic 
as developed in this report will call for a many turn spiral 
approach. 


1 . 4 Importance of Analytic Approach 

This system has as its central feature a turn 
capability prediction based on a simple analytic 
formula. That is, based on the present state of the 
vehicle and an estimate of the vehicle lif t-tordrag 
ratio, L/D, a turn capability of say 2 1/2 revolutions 
is predicted at the start of circling flight. The guidance 
will then choose a commanded bank angle based on the nominal 
value such that desired number of, say two, revolutions 
result; and this computation is repeated every cycle. 
Subsequent analysis will show quite close agreement between 
the actual and predicted turn capability. 

This analytic type approach is important for several 
reasons. First, a natural type trajectory results. That 
is, the vehicle is not commanded into any extreme maneuvers 
by the very fact that the guidance is based on its predicted 
capability. Second, the simple nature of this approach 
allows a flexibility to adapt to changing constraints that 
always appear as the program matures. This is particularly 
important in the early stages of the program. Third, the 
simple nature of the program accelerates the check-out 
and debugging time and also does not over-tax the airborne 
computer storage or duty cycle. This last feature is of 
course quite secondary in view of the vast computer capability 
proposed for the Space Shuttle. 

It will be demonstrated in section 4 the analytic 
techniques are possible for both the transition and flare 
maneuvers as well as the energy dissipating turn. Further, 
it will be shown that several approach geometries are 
possible with this analytic technique even though it was 
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primarily designed for circling flight. 


1.5 Important Conclusions 

The most important lesson that will be demonstrated 
in subsequent chapters is that simple analytic formulae 
will meet all mission requirements. Moreover, this type 
approach is highly desirable since new requirements can 
be easily adapted by the analytic nature as they may 
arise. The guidance permits ranging to all points 
within the vehicle capability. Since important 
disturbances such as winds are compensated directly, 
there is a high tolerance to disturbances. 

Key vehicle parameters such - as L/D need not be 
estimated in flight except perhaps for the final flare 
maneuver. Since this estimation is an easy task, there 
is little reason for not doing it. 

The large initial navigation errors are resolved 
within the first few measurements and there is weak 
coupling between the navigation errors and the steering. 
However, large cross-track errors can lead to navigation 
filter divergence if measurements are not available to 
resolve these type difficulties. That is, there is a 
real danger of underspecifying the navigation equipment 
to the point of failing the mission. 


Preceding page blank 
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CHAPTER 2 

VEHICLE CHARACTERISTICS AND CAPABILITY 


2.1 The 040a Vehicle 

The Shuttle vehicle for this study is the 040a 
model which was the baseline model for many MSC 
studies in this period. The aerodynamic characteristics 
are defined in ref .ROl . The computer model based on 
this reference is found in Appendix D. Typical L/D values 
are shown in Figure 2.1 and more detailed characteristics 
(C L C D & L/D) for the subsonic flight are shown in 
Appendix F. 

It is seen that the "clean" vehicle has a maximum 
L/D of 7.6 occurring at an angle-of-attack of 8 degrees. 
Note also that the speed brake has a profound effect 
decreasing the L/D almost in half when fully deflected. 

The nominal vehicle has these other significant 
characteristics : 

WT 1400001b. 

Ref Area 3130 ft 2 

Mean Chord 51 ft. 

The ground effect from ref. El o3 is also simulated. 
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Fig. 2.1: 04 0a L/D Characteristics at Mach .5 


2.2 Footprints 


The subsonic ranging capability has been defined 
in Ref. [1] . This capability is shown in the form of 
the- footprints in Fig. 2.2. Here are shown range capability 
for a vehicle with L/D=8, starting at various initial 
altitudes. The negative down range direction is smaller 
because two 180 degree turns are required to attain this 
point. The first changes the heading to the negative 
direction. The second changes the heading back to aline 
with the runway direction. It is seen that near 80 nm 
are possible in the forward direction, with 50 nm in the 
negative direction if one starts at a nominal 60,000 ft. 

The supersonic ranging capability is considerably 
less because the maximum L/D is only 2. That is to say 
variations in range are less pronounced because there is 
less possible variation in L/D. It is o£ interest that 
the glide slope is not even close to 1/ (g-) at this low 
L/D. The altitude rate can be foundry using the results 
of Appendix A by first noting that V z g is much greater than 

2g 

1 in Eq. (A. 13) so that assuming <j> = 0, we have 



Then Eq. (A. 4) gives for dH 

dt 


( 2 . 1 ) 


dH 

dt 


2g 

<£>ve 


( 2 . 2 ) 


Altitude rate thus varies from -120 fps when V = 6000 fps 
to -720 fps when V = 1000 fps. 


The constant deceleration Eq-. (2.1) of about 1/2 g is 
appropriate to generate the supersonic range as has been done 
in the transition ranging logic (Section 3.2, Eq. 3.2). 
Applying Eq. (2.1) into Eq. (3.2), gives a maximum range of 
180 nm for the supersonic phase, assuming an — of 2. 

No supersonic footprints have been generated here. 

This footprint capability is rather shown in the performance 
data of Section 4.3, where the guidance achieves edges of the 
footprint by flying at maximum L/D. 
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CROSS RANGE ( N Mi ) 



Fig. 2.2: Subsonic Range Capability of Shuttle Orbiter (L/D=8) 




CHAPTER 3 


GUIDANCE DESIGN 


3.1 Guidance Phases 


The guidance designed here was essentially completed 
as part of the Phase B effort [1 ]. Several improvements 
and additions were made since that time including transition 
logic, improved targeting, damping terms added and flare 
logic. Also, a major effort was made in testing the guidance 
in a realistic navigation environment. For sake of completeness, 
the total guidance logic is described here. 

The trajectory starts at about 150,000 feet with the 
vehicle at Mach 6, having, just emerged from radio black- 
out. There is then a transition maneuver from high angle- 
of-attack of 30 degrees to a low value of about 8 degrees. 

During this ' transition maneuver, there is some ranging also 
performed though only a limited amount is possible. Next 
there is a glide to the pre aim-point before the runway, 
performed at maximum L/D. Circling flight is next 
initiated with one or two revolutions possible. The next 
phase is the final approach starting at 10,000 feet 
altitude and ten miles before the runway threshold. A 
first flare maneuver starts at roughly 700 feet altitude 
to reduce the glide slope from 10 to 3 degrees. Finally, 
a second flare is initiated at about 50 feet altitude to 
reduce to sink rate to a nominal 4 feet/sec. 


3 . 2 Transition Logic 

There are two functions accomplished during the 
transition phase: 1) a change in angle-of-attack , 

2) ranging to an aim point near the landing site. 

The change in angle-of-attack can be accomplished 
with a pitch over rate of between .05 and .2 deg/sec 
without either exceeding the g limits on one hand and not 
accomplishing the maneuver in time on the other. There is 
therefore, established a velocity at which to start the 
maneuver, designated VLM2. There after, the angle-of-attack 
is commanded as a function of velocity by: 

a = a o - KTR (VLM2-VA) (3.1) 


where 


a = angle-of-attack 
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a = initial value 
o 

KTR = Pitch over gain nominally .0075 
VA = velocity relative to the ground 


The ranging portion of the transition maneuver is 
targeted to a point 30 nm before the runway with an arrival 
velocity of 1000 ft/sec. This point is roughly centered 
within the remaining ranging capability of the vehicle. 

The roll angle is commanded to accomplish this ranging 
by first calculating a nominal range based on an assumed 
constant rate of change of velocity. This constant 
deceleration has been observed over a large class of 
transition maneuvers and leads to a simple range calculation: 

RTOGON = (VA 2 - VF 2 )/2ACTR (3.2) 


where 

nominal range 
final velocity = 1000 fps 
nominal deceleration 
15 ft/sec 2 

The roll angle is then commanded by 

<P = cos' 1 (RTOGO/RTOGON) (3.3) 


RTOGON = 
VF = 
ACTR = 


where 


RTOGO = range to initial target from Eq. (3.6) 

$ = commanded magnitude of roll angle 

The roll command is modified with a lateral logic which 
will allow only one roll reversal and turn-rate limiter which 
limits the azimuth rate so as not to aggravate the 
effect on the sonic over pressure (“sonic boom"). The roll 
angle is initially commanded in the direction so as to 
decrease the heading error . (Or it would be initially 
commanded in the same direction as in the previous entry 
guidance phase) . This roll direction is maintained until 
the heading error, DEL, exceeds a threshold value, DELTR. 

Then, the opposite direction of roll is commanded. 

The turn rate limiting is not started until velocity 
is less than 3000 fps, (VLM3 ) . Excessive turn rates are 
not possible at higher velocities. At this time, the turn 
rate at nominal roll angle is calculated by: 

n = g tan(|^) (3.4) 
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where 


fi = turn rate 
rlb = nominal bank angle 
g =32.2 ft/sec 2 


If the turn rate exceeds the limit value, OMLM, the roll = 
is commanded to achieve this limited value by: 




^ -1 ,OMLM VA, 

tan ( ) 


(3.5) 


Also, the roll angle is commanded to drive the heading angle 
to zero in the next interval if this is possible and no 
limits are exceeded. 


3.3 Landing Pattern Targeting 

When the velocity is less than VLMl (1200 fps) , the 
start of the circling flight mode is initiated. An aim- 
point is selected to start the turning flight so as to 
arrive at the final approach with no lateral offset by 
the following two step procedure: First, the target vector, 

URT, located at the edge of the landing field is projected 
backward a fixed bias distance, to UR'T2. Figure 3—1 is 
appropriate. Ten miles has been selected for this bias 
distance in the initial studies. This bias corresponds to 
an altitude of 10,000 feet at the start of the final glide 
and a nominal L/D of 5 for the final glide. Then using the 
heading of the vehicle to this initial aim point, an offset 
aim point is calculated, see Figure 3-2. Two circles are 
calculated. The circle of radius R corresponds to circling 
flight at the velocity V! which is the arrival velocity 
at URT 2 . The circle of radius r corresponds to circling 
flight at V 2 , the equilibrium velocity at the final altitude, 

H f . The calculation_f ollows . First calculate the range 
from the vehicle to URT2: 

RTOGO = cos" 1 (UR.URT2) SF ■ (3.6) 

where 

RTOGO = range to offset target, feet 

UR = unit vector at present position 

URT2 = initial target unit vector 

SF = scale factor converting earth central angle 
to feet 

Then calculating the altitude loss 6H in this segment assuming 
flight at a known maximum L/D, (L/D)^ 

6H = RTOGO/ (L/D) m (3.7) 
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VEHICLE 



Fig. 3-1 Initial Target Calculation 


_ VEHICLE UR VECTOfi DIAGRAM REFERRING TO VELOCITY 



Fig. 3 -2 Final Target Calculation 
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The velocity at URT2, Vi, is then assuming equilibrium at 
both the initial and final points. 

V x = VA1 e- (P 6H/2> (3.8) 

where 

VAl = present velocity corrected for roll angle 
= VA V g/AMAG cos (<j>) 

g ss reciprocal of atmospheric scale height 
= 1/25000 

AMAG = total acceleration, fpss 

Similarly the velocity at the start of the final glide, V 2 , 

is related to the altitude at the start of final glide, H f , 

by 

V 2 = VAl e- 8 ( h -H£)/2 (3.9) 

The initial and final turn radii (R,r) are then 
calculated using the commanded bank angle <f>. 

R = Vi 2 /(g tan <j> ) (3.10) 

r = V 2 2 /(g tan ^ (3.11) 

where 2 

g = gravitational acceleration = 32.2 ft/sec 

= cf>n initially, then defined by Eq.(3.22) 

4> = nominal bank angle, 30 deg. 

The final offset vector is then seen to be 

URT3 = URT2 + ( r " R R g° S ■■ — ) URTl + R 9 UF (3.12) 

where 

URT3 = final offset unit target vector 
URTl — UF x URT 

URT = unit target vector at edge of runway 
UF = unit vector in direction of runway 

The radius of the earth, RE, is used to convert the offset ' 
distance into angles. The angle 0 is not calculated explicitly, 
but the cosine and sine of 0 are calculated by the vector 
relations 

cos 0 = UN I1 > URTl (3.13) 

sin 6 = -V 1 - cos 2 6 SIGN (UNIl • UF ) (3.14) 


where 
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UNI1 = UNIT (URT2 x UR), for initial calculation; 

“ UNIT(URT3 x UR) , subsequently 

The appropriate vectors are shown in Figure 3-1 and SIGN 
is the sign function. 

The above calculations are used for the initial 
calculation of the offset target vector', URT3. A precise 
calculation would involve an iteration since the target 
vector is a function of the range-to-go, which is in turn 
a function of the target vector. This iteration is avoided 
by repeating the calculation of the offset target vector 
on subsequent guidance cycles replacing URT2 with URT3 in 
Eq. (3.6). The offset target calculation stops when the 
heading error, Del, is less than .01 radians. 


3 . 4 Initial Turn 

The initial turn is accomplished by banking at the 
nominal bank angle to drive the heading error toward zero 
until a calculated reduced bank angle will drive the error 
to zero in the next time interval . 

The heading error Del is 

Del = cos -1 (UNI1 • UNI) (3.15) 

where 

, UNI = UNIT (V x UR) 

The bank angle command is 

4> = 4>n SIGN (UNI3 • UNI)* (3.16) 

where 

UN I 3 = UR x UN II 

If ftdt is less than DEL, the turning rate ft is calculated 
by 

ft = g tan <}>n/V (3.17) 

and dt is the sampling interval. Five seconds are used in 
the simulations. 

If ftdt is greater than Del, the bank angle is commanded by 
<t> =' tan^ 1 t ^ - 1 - ^ ) SIGN (UNI3 ♦ UNI) (3.18) 

* Note that positive bank directs a component of lift along UNI. 
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Throughout the initial heading change the angle- of - attack 
is kept fixed near the value for maximum L/D. 

The above steering logic neglects the dynamic response 
of the vehicle to the roll commands and no dynamic lag has 
been included in the simulations thus far. It is expected 
that the "dead beat" response described above can be extended 
to include these dynamics and would be superior to a constant 
gain feedback system type design. This is because more 
information is available on the vehicle characteristics 
and the implementation is more suited to the digital 
computer needed for guidance and navigation calculations. 


3 . 5 Glide to First Aim Point 

The glide to the pre-aim point is accomplished at 
maximum L/D. During this interval a calculation is made 
of the turn to lose altitude required during Phase 3. This 
turn angle is 

D0 = (L/D) sin <J>n g HS ^-)RTD (3.20) 

v 2 v x 

+■ (L/D) m sin 4>n (H-H f ) RTD/2HS 

This result is derived in Appendix A. Somewhere between 
1 and 2 revolutions will be required depending on the range 
to the pre-aim point. 

The next phase is entered when the vehicle^ arrives at 
URT3. This is detected when the dot product, UNIl • UNI, 
is negative. 


3. 6 Circling Descent 


3.6.1 Turn Angle Control 

The nominal circling descent is at constant bank 
angle <i>n. This angle has been set to 30 degrees for 
initial studies. The turn angle is calculated using 
Eq. (3.2 0) with V used in place of V1-. The required turn 
angle is then calculated as 

D61 = n 360 + Daz (3.21) 


where 

Daz = Az - Azr 

Az = Azimuth = cos (UE • UNI)RTD 

if V*UE is negative, Az = 360-Az 
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Azr = Azimuth of runway 
n, chosen so that |D0 - D0l|<1.8O 
RTD = factor to convert radians to degrees 
UE = unit vector in easterly direction 

An option to the logic allows n to be an input variable 
if desired. With this option, one can select the desired 
number of turns but only, of course, before the actual turn 
has started. 


The number of revolutions, n, is further modified so 
the value of D01 does not differ by more than the 18 0 
degrees from that of the previous computation cycle . . 

Without this feature, an induced oscillation is possible 
as the desired turn angle changes back and forth to values 
near - 180 degrees from the predicted turn angle. The 
bank angle is then 

<}> = sin ^ (3.22) 

and the bank angle is limited to values less than <f> max r 
now chosen at 45°. 

Thus, we see that the bank angle is continually 
adjusted so that the heading is the correct value at the 
final altitude. The initial targeting assured that the 
position relative to the runway is correct when this final 
heading is achieved. 

A further modification of the turn logic allows for 
straight flight segment when the desired turn angle is 180 
degrees. This logic is normally disabled because improved 
initial targeting does not require it. But it is included 
here to allow trajectory shaping as has been done in Chapter 4. 
This is done so that the predicted turn angle can decrease 
to the desired value. For example, the situation could arise 
where the desired turn angle was 180 degrees and the 
predicted turn angle, 360°. In this case the commanded 
bank angle would be 14.5 degrees and a much larger than 
nominal turn radius would result. The altitude loss during 
this straight flight segment is added to the final altitude, 

Hf . This is done to allow for the additional range needed 
in the final phase. Typically, one or two nautical miles _ 
will be flown in the straight-flight portion and correspondingly 
about 1,000 feet will be added to the final approach altitude, H f . 
In this situation, however, the improved targeting would allow 
a 14.5 bank angle turn and no straight segment would be required. 
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3.6.2 Trajectory Damping 


In the transition, the turn and the glide to the 
initial aim point, the phugoid motion is damped by 
commanding an incremental angle-of-attack. 

Da = KDMP (Rref-R) 

where 

Da = incremental angle-of-attack 
KDMP = constant damping gain 
R = altitude rate 
Rref = reference altitude rate 
= - 250 fps in transition 
= -VA/{L/D) m , before turn 
= -VA/((L/D) m cos <|>) , during turn 

The angle-of-attack is further modified to 

anticipate the transient caused by changing angle-of-attack 

commands by . 

{cos <j>) Q 


where 


(cos <j>) = cos <j> from last computation cycle 


3.7 Center of Curvature Control 


To make a positive control of position relative to 
the landing field during the turn, logic is written to 
control a specified center of curvature. This logic 
would be appropriate for the case of disturbances such 
as wind, improper navigation and unknown vehicle variations 
since the initial targeting assured proper positioning in the 
no disturbance case. The logic adds an incremental roll angle 
to the commanded roll angle so as to decrease distance 
between the instantaneous and desired center of curvature. 
Referring to Figure 3-3 a unit vector at the desired center 
of curvature, Uc, is described by 

Uc = UNIT (URT2 + ^|- URT1) (3.23) 

where r is defined by eg. (3.6). 

The instantaneous center of curvature, Ucl, assumes turn 
at nominal bank angle 

_ DT 

Ucl = UNIT (UR + ~ UNI) (3.24) 


where 


R1 = V 2 / (g tan 4>) 
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The distance, Dele, between the desired and instantaneous 
centers of curvature is then 

Dele = cos" 1 (Uc • Ucl) SF (3.25) 

The incremental roll angle, DRL, driving Ucl toward Uc 
is commanded by 

DRL = — KCN Dele SIGN (CTEST) (3.26) 


KCN = constant gain = .0003 
CTEST = UNI • Uc2 

Uc2 = UNIT (Uc x Ucl) 

= normal to plane containing Uc and Ucl 

Note that this incremental roll^ angle has the_effect of 
increasing the turn radius if Uc is ahead of Ucl^ (As in 
Figure 3-3), and decreasing the turn radius of Uc is behind 
Ucl. In both cases, the instantaneous center is driven 
toward the desired center. 


3. 8 Final Glide 


The final glide starts when the indicated altitude 
is less than Hf . At this point the target is repositioned 
at the edge of the runway, at URT. Heading corrections 
are then made with logic identical to Eqs. .(3.10, 3.11, 
3.12) replacing only URT3 with URT. 

An "s-turn” logic is exercied so that the lateral 
position as well as the heading is controlled. The lateral 
displacement accrued, d, in correcting the heading error 

is calculated by 

« 

d = 2 Rl sin 2 (— ^— ) SIGN (7LAT ) (3.27) 

where 

Rl = turn radius at nominal bank angle 
= V 2 /(g tan $n) 

d = lateral displacement during turn 
(positive if along URTl) 

VLAT = V • URTl 

The lateral error, Dl, is calculated by 

D1 = (UR ♦ URTl ) SF (3.28) 

If the total lateral error, Dl + d, is greater than a 
threshold value, the heading error is modified by 

Del = Del - Dei 3 SIGN(d) (3.29) 


where 


Del- =45 degrees nominally, but is varied by Eq. 

J (3.32) to account for long range cases. 
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URT2 


URT . 



Fig. 3 -3 Geometry for Center of Curvature Control 


URT 



Fig. 3 -4 Geometry for "s-turn" Logic 
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• This logic has the effect of heading so as to close the 
lateral error and then reversing the bank angle to correct 
the heading error at a time when the lateral error .will 
be zero when the heading is correct. Only one "s turn" 
is allowed in final glide. A standard linear feedback 
control on heading error and lateral displacement is then 
initiated after the initial "s-turn.” 

The. .law. commands—the- bank-angle -by- terms— proportiona-1 — 

to the lateral displacement and lateral velocity 

<j> = KDD D1 + KW VLAT (3.30) 

where 

KDD = -.07 = position gain 
KW - -.36 = damping gain 

The gains have been chosen to provide a closed loop natural 
frequency of one minute and a damping ratio of .7. 

At this point the angle-of-attack is modulated so 
that flight path intersects URT. This logic will be made 
with modulation on the "front-side" of the L/D curve. The 
angle of attack is commanded by the following equation: 

K (L/D - L/D ) 

“ = “o + 5 (L/D) /Oa + KDMP 0*ef-R> (3.31) 

where 

a Q = 6 degrees 
L/D o = 5 

9 (L/D) 3a =1.4 deg -1 
K = 2 

a is limited between 4 and 15 degrees 
L/D = (RTOGO-BIASl 6080) /H 

C 

Biasl = Bias on final phase range = 1.8 n.m. 

Rref = Reference altitude rate = -VA/{L/D). 

KDMP = damping gain = . 01 

3.8.1 Variable Lead Angle in S-Turn 

Should the range to the landing site be so long 
that no turn is possible, special logic is in order. 

Because the conditions can arise that one turn is too much 

and a striaght-in approach is too little, an overshoot 

may result. In this case, the final phase is entered directly 
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if the predicted turn angle is less than 180 degrees, and 
the excess energy between one half turn and straight—in is 
dissipated by varying the lead angle DTHl in the s turn logic, 

L/D - L/D 

Del 3 = Del 3Q (1 ~~~2 (3.32) 


where 

Del 3Q . — .75 rad. 


3.8.2 Speed Brake Law 

The Speed brake is needed mainly to decrease the 
reference L/D during the final approach while still retaining 
a reasonable lift coefficient. Refering to Fig. 2-1/ we see 
that L/D of (the reference value) 5 for the clean vehicle is 
at a 3 degree angle-of -attack. This low angle-of-attack 
would result in a high velocity to maintain equilibrium because 
the lift coefficient is so low, while a speed brake deflection 
of 25% doubles the anlge-of -attack for the reference L/D of 5. 

The speed brake is used primarily to control velocity 
by defining a reference velocity as a function of altitude 

VRSB = VR + KSBH H (3.33) 

o 


where 


VR = reference velocity at zero altitude 
o 

= 425 fps 

KSBH = speed brake gain = .003 

The speed brake is then deflected to achieve this 
velocity by 


SB = KSB (VA-VRSB) + SBO (3.34) 

where 

SB = speed deflection in per cent 
KSB = speed brake gain = .5 

SBO = bias setting = 20% 

A limited amount of testing of an additional term to range with 
the speed brake showed little improvement and was discarded. 
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3 . 9 Flare Control 


The flare maneuver is compo'sed of two segments. 

The first flare reduces the flight path angle from about 
10 to about 3 degrees. The second reduces the sink rate 
to a nominal 4 ft/sec. This logic commands a constant 
incremental normal acceleration to accomplish both flare 
maneuvers. The first flare is with an incremental g of 
nominally 0.3 g initiated so as to be completed at an 
altitutde of HF2 normally 600 feet. To do this, a 
threshold altitude is calculated at which to start the 
maneuver so as to be complete at HF2. 

HTT = HF2 + (R 2 - *^)/2 GF2 (3.35) 

where 

HTT = threshold altitude to start maneuver 
R 2 = final altitude rate = 18 fps 
GF2 = incremental g = 10 fpss 
HF2 = altitude at end of flare 

When the altitude is less than HTT an incremental 
g is commanded by increasing the angle-of -attack 

a = a (GF1 + GS)/AMAG (3.36) 


where 

• « 

GFl = (— R + R 2 )/dt, limited to GF2 

When H is less than the flare altitude, HF2, ranging is 
accomplished by commanding the altitude rate by 

R 2 = - VA H/ (RTOGO - BIASFL) (3.37) 

where 


BIASFL = aim point before runway 
= 2200 feet 

Then angle-of -attack is commanded by Eq. (3.36). 

The second flare is started when the altitude is 
less than the threshold value (HTT2 =50 ft) and is accomplished 
by commanding angle-of-attack with Eq. (3.36) replacing the 
limit incremental g with GF3 = 6 fpss. 
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3.10 Left Hand Turn Logic 


Originally, left hand turns were desired to possibly 
handle the long-range cases where no full turn was possible. 
But these possibilities were handled by the variable lead- . 0 
angle logic of Section 3.7.1. However, other constraints 
may require this type turn so the modifications to the above 
logic are described here. There are four. 

First, re-target the offset vector by modifying 
Eq. (3.12) 

URT3 = URT2 - ~ r E C ° S ^ ^RTl + R 6 UfJ 

^ (3.125.) 

Second, change the sign of the azimuth error 
Daz = - (A z - A zr ) (3.125.) 

Third, command the opposite signs bank angle 
<J> = - sin -1 (^j-) sin <j>n) (3.225.) 

Finally, redefine the unit vector Uc by 
Uc = UNIT (URT2 - URT1) (3.235.) 


3.11 Concluding Remarks 


The above guidance scheme can be found twice more 
in this report in increasing levels of detail. First, 
detailed flow graphs are presented in Appendix B. Second, 
a Fortran listing of this program is found in Appendix D. 
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CHAPTER 4 

GUIDANCE PERFORMANCE 


4.1 Nominal Case 


The logic was tested with a variety of initial 
conditions and landing sites. The nominal trajectory 
is described first. This trajectory starts near 
150,000 feet just prior to the transition maneuver. 
The landing site is chosen to be centered within the 
vehicle capability from this point. The vehicle 
initial state is summarized in Table 4-1. 


V 

= 

6055 fps 

H 

- 

143020 ft. 

Y 

=z 

-1 . 9 deg 

azimuth 

= 

45° 

Mach 


5.71 

a 

- 

29° 

latitude 

= 

0 

longitude 


0 

Runway latitude 


2.4 deg 

Runway longitude 

= 

2 . 4 deg 

Runway azimuth 

= 

75 deg 


Table 4-1: Vehicle Initial State and Nominal Runway 


The ground track for this trajctory is shown in 
Fig. 4.1. On this trajectory are shown tick marks for 
each minute of flight. Almost 17 1/2 minutes elapse 
(1040 sec) . Shown also are the initial pre-aim, and the 
target for transition maneuver, as described in Section 3.2. 
The transition maneuver starts at about 9 0 seconds -and ends 
at the six minute point. In this trajectory, a one-turn 
circling flight before the runway is performed. 
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An expanded view of the ground track and an 
elevation view of the last 10 minutes of flight 
are shown in Fig. 4.2. Shown here is the offset 
target and desired center of curvature. 

In Fig. 4.3 are the altitude, velocity and 
altitude rate time histories. Altitude rate is seen 
to be about -250 feet/sec during the transition 
maneuver and shows small peaks at the start of both 
the turn and the final approach. Steps in altitude 
rate are also seen for each of the flare maneuvers. 

The altitude here shows a 2 Slope character: Steeper 

for the first 400 seconds due to the supersonic L/D 
or 2, and later shallower with the subsonic L/D of 7. 

The velocity decreases in a monotone fashion with the 
sharp decrease at 1000 seconds accompanying the flare 
maneuver . 

The control variables in Fig. 4 .4 show the roll 
angle, angle-of-attack and speed brake deflections. 

The roll angle shows the reversal at 150 seconds as 
part of the transition ranging. At 360 seconds there 
is an impulse of roll angle as the vehicle heads toward 
target 2. Circling flight begins at about 600 seconds 
with a near constant bank angle of 20 degrees. At 900 
seconds the turn is complete and the final approach begins. 

The angle-of-attack starts at 26 degrees (not 29) 
due to the phugoid damping which attempts to maintain an 
altitude rate of -250 ft/sec. At 90 seconds the pitch over 
is begun, though this is partially masked by the damper 
actions. The transition is complete at about 400 seconds 
arriving at an angle— of-— attack of 7 degrees. The final 
approach shows an alpha of 6 1/2 degrees starting at about 
900 seconds. Finally, the flare maneuvers show the double 
pulse character for the first and second flares and charac- 
teristic rise needed to maintain “equilibrium" as the vehicle 
slows down. 

The speed brake shows a pulse to about 40% at the start 
of the final approach finally settling to a 20% deflection. 

Not shown is the landing gear deflection which 
occurs at the start of the final approach. 
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The dynamic pressure in Fig. 4.5 is not constant 
during the turn as the turn angle prediction implicitly 
assumes. This is because of the model of the 040a vehicle 
which shows a slow decrease of C L with Mach number. Dynamic 
pressure must therefore increase to maintain equilibrium. 

The characteric drop as seen in dynamic pressure at the end 
as the vehicle is pitched up on the flare maneuver. , 

The total acceleration, also in Fig. 4.5, does not 
stray far from the one g level. The peak is 1.25 g. The 
double pulse corresponding to the two flares is also seen in 
the acceleration trace as it was in the angle-of-attack trace. 

The predicted and desired turn angles are shown in 
Fig. 4.6. It is seen that these two angles converge as 
driven by the 20 degree bank angle. A nonlinearity has been 
noted in the turn angle to bank angle. Namely, increasing 
the bank angle from the nominal 30 degree bank angle shows 
less of a change in turn angle than does decreasing the 
bank angle. 

Details of the flare maneuver are seen in Fig. 4.6 
showing altitude, altitude rate, velocity and angle-of- 
attack as a function of range-to-go to the landing site. 

The first flare shows an step in alpha of about 2 degrees 
lasting about 2 seconds. The second flare lasts only 
about 3 seconds. There is a characteristic drop in alpha 
just before touchdown as the ground-effect becomes effective. 

The touch down velocity of 261 fps is low. This is 
made possible by using the ground effect which allows an early 
flare. This should probably be increased by delaying the first 
flare (decreasing HTT) and at the same time decreasing the 
bias on the final phase (BIAS1) . 
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Fig. 4.1: Ground Track for Nominal Case 
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Fig. 4.5 : Acceleration and Dynamic 


Pressure Time Histories for Nominal Case 
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Fig. 4.f: Altitude, Altitude Rate, Veil 

as Function of Range to go D 


Nominal Case 


4.2 Variations on Nominal Case 


4.2.1 Two Turn Case 

A two turn case for the nominal target is shown 
in Figs. 4.8 4.11. This is the more natural mode for 

the nominal target centered in the footprint. It was 
generated by setting the angle-of -attack after transition 
to the max L/D value of 8 degrees and letting the steering 
equations choose the desired number of turns. Only the 
final part of the trajectory starting at 400 seconds is 
shown . 


The ground track is in Fig. 4.8, There is a 
characteristic drift toward the landing site from the 
desired turn center. This is due in part to the 
increasing roll angle during the turn, but more directly 
it is caused by the fact that a higher velocity and slower 
turn rate at the start of the turn would naturally 
contribute to such a drift. 

Altitude, altitude rate and velocity traces 
in Fig. 4.9 are very similar to the nominal case 
transients marking the beginnings and end of the turn 
segment are seen in altitude rate. Small transients at 
700 and 750 seconds also mark the step changes in roll as 
the center control is exercised. 

The control variables in Fig. 4.10 are also similar 
to the nominal case. Angle-of -attack is near constant near 
8 degrees until the start of the final phase. 

Roll angle starts near the nominal value of 30 degrees, 
but increases slowly as the initial turn angle predictions 
overpredicted the turn capability, a step down, then up, then 
down is seen in the roll angle as controlled by the center 
of curvature control logic. 

The speed brake shows a nearly constant 30% deflection 
starting at 850 seconds, the beginning of the final phase. 

The dynamic pressure in Fig. 4.11 again shows its 
characteristic rise, this time partly in response to the ever- 
increasing roll angle. The acceleration does not wander 
far from the 1 g level. 
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4.2.2 "Race Track" 

In response to the discussion with pilots [ 5 ] 
and the results of other MSC investigations [ 4 ] f a two 
turn trajectory was generated with this logic by changing 
certain control parameters. This would appear to be a 
back-up mode for this guidance should the pilot have to 
take-over from a primary system failure. The roll angle 
was constrained to 30 degrees and the target point was set 
at 5 nm before the runway and the corresponding final 
approach altitude set to 5000 feet. In this way, the vehicle 
would pass over the runway on the first turn. 

The resulting trajectory is shown in Figs. (4.12) 
through (4.15). The ground track in Fig. (4.12) shows 
the characteristic down wind leg away from the runway. Though 
the runway is only about 20 nm away at most, this may not 
be desirable for the nominal landing. Note that the initial 
targeting provided an initial turn radius of the right value. 

The altitude rate in Fig.. (4.13) is "busier" in this 
trajectory as turns start and end. But this is not particularly 
ob j ec tionable . 

The control variables shown in Fig. (4,14) show 
slightly more variation in angle-of -attack. The constant 
roll angle segments might be easier for a nominal mode. 

The acceleration and dynamic pressure in Fig. 4.15 
also show more variation than do the other two cases. 


4.2.3 Left Turn Case 

A left turn nominal case is illustrated in Fig. 4.15a. 
This is in response to the guidance changes described in 
Section 3.9. It is expected that right hand turns would be 
used universally. This would make the pilot monitor task 
easier. Only would an unusual mission constraint, not yet 
defined, call for a left hand turn. 
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Fig. 4.15a: Ground Track for Nominal Target with Left Turn Logic 



4*3 Different Landing Sites 


Landing sites other than nominal were tested to see 
the capability of the guidance id range to the capability . . 
limits of the vehicle. This of course is the same as varying 
the initial position of the vehicle relative to a fixed 
landing site. Such variations might occur due to faulty nav- 
igation during the preceding trajectory segments. A once 
around. abort mission with no navigation update has been seen 
to give particularly large initial errors. In this case/ 
errors on the order of 30 nm are possible. 

Ten different landing sites are shown in this section. 

In all cases, the vehicle starts in the same state as 
described in Table 4.1. The landing sites are in the order 
of decreasing range with two large cross range sites 
listed last. 



■Landing Site Location 


Case 

Latitude 

(deg) 

Longitude 

(deg) 

Comment 

1 

3.3 

3.3 

- 

2 

2.6 

2.6 

- 

3 

2.4 

2.4 

- 

4 

2.2 

2.2 

- 

5 

2.0 

2.0 

- 

6 

1.5 

1.5 

- 

7 

1.0 

1.0 

azimuth rate 
limit off 

8 

.5 , 

.5 

missed 

9 

1.5 

3.5 

large crossway 

10 

3.5 

1.5 

large crossway 


Table 4.2 Landing Site Locations 
Runway azimuth = 75° 
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In most cases, a two turn trajectory results because 
at this point a one turn trajectory had not been selected 
as nominal and a two turn trajectory was within the vehicle 
capability. 

Case one, target at (3. 3, 3. 3), shown in Fig. 4.16, 
represents the extreme long range case. An early transi- 
tion starting point was ev.en required to squeeze out the 
extra 10 nm that this allows. 

The circles increase slightly in diameter as the 
range to target decreases from case to case until we 
come to the target at (1.5, 1.5). In this case, the 
guidance must try hard to get back to the final approach 
path and only a half turn on to final approach is possible. 

The' next case, target at (1. 0,1.0) Fig. 4.22, is not 
possible with the azimuth limiter on. For that matter, it 
is near the minimum possible range. An approach from the 
other runway directions should be considered in this case. 

The effect of the variable lead angle logic. Section 3.8.1, 
is seen starting at 1.5 degrees latitude. 

The landing site at (.5,. 5), Fig. 4.23, is not attainable 
with this logic. It appears to be outside the vehicle 
capability. 

The two large cross track landing sites are shown in 
Figs. 4.24 and 4.25. They are ± 85 nm from initial 
trajectory plane. Somewhat greater lateral range is 
possible but was not tested. 

In all, landings are demonstrated within a band 
195 nm in the down range direction and 170 nm in the 
cross range direction. 
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Figure 4.16 Ground Track Target at (3.3, 3.3) 
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Figure 4.17s Ground Track Target at {2.6, 2.6) 
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Figure 4.19: Ground Track for Target at (2.2, 2.2) 
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Figure 4.20: Ground Track Target at (2.0, 2.0) 
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Fig. 4 * 21 : Ground Track for target at (1.5, 
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4.4 Wind Performance 


4.4.1 Semi Steady Winds 

That the design wind for Shuttle [ 15] is an 
extreme test, is shown in Fig. 4.26. Here the Shuttle 
attempts to land into a head wind with no compensations.. 

The profile of this wind is also shown in Fig. 4.26. 

The guidance described in Section 3 cannot 
compensate for this wind as shown. Successful landings 
with winds of 1/2 this magnitude have been made, however. 

The center-of-curvature control is the reason for this 
capability. The vehicle is blown off course and the center- 
of-curvature control brings it back. 

Three additions to the logic will allow this capability 
however. They all involve knowledge of the vehicle air speed. 
Heretofore we steer using fround speed as inferred from the 
inertial measurement unit. 

The first is to replace VA with VAIR, the air speed 
in the steering equations. Second change the bias on the 
landing field at the start of the offset target calculations 
by 


BIAS = BIAS + KBIAST (VA-VAIR) 


(4.1) 


where 


KBIAST = .03 nm/fps 

Note that a head wind will decrease the bias while a tail wind 
will increase it. 

Third, modify the commanded angle-of -attack in the 
final phase by using this difference by an amount Da. 

Da = KWF (VA-VAIR) (4.2) 


where 


KWF = . 01 deg/fps 


With these modifications, the trajectory for a 
head wind and cross wind result as shown in Figs. 4.28 and 
4.35. Of most interest is the extreme action of the roll 
angle in Fig. 4.30 and 4.34, in compensating for the wind 
action. Note also that air speed not ground speed is 
displayed in the velocity plots. 
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Fig4,2'6i: Ground Track & Elevation Views for CE3cfcre*ae Head Wind 
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4.4.2 Gusts 


Only limited testing with gusts was made. A gust 
of the 1-cosine form was introduced during the final 
approach phase where the gust effect is most pronounced. 

These tests are illustrated for a downward gust 
and a cross wind gust introduced at 950 secs and lasting for 
30 seconds. The peak value of the gust is 40 ft/sec. The 
down wind gust effect on the nominal trajectory is 
in Fig. 4.35 and the cross wind gust is in Fig. 4.36. 
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Fig. 4.'>2B: Altitude Velocity & Altitude Rate Time Histories for Headwind Case 
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Figure 4. 29: Roll, Alpha, 
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Fig. 4 .30: Acceleration and Dynamic Pressure Time Histories for Headwind Case 
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Fig* 4*34: Acceleration and Dynamic Pressure Time Histories for Cross Wind Case 
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Fig. 4.35 : Ground Track for Down Gust During Final Approach 
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Fig .4. 3 6 : Ground Track for Cross Wind Gust During Final Approach 



4.5 Vehicle Parameter Variations 


In this section the vehicle mass and the vehicle 
L/D are each varied ± 20% for the nominal trajectory. 
Satisfactory performance is seen in all cases except 
that more information is needed in the flare logic in 
the low L/D case. 

The heavier vehicle ground track is shown in 
Fig. 4.37. A higher velocity corresponding to the 
higher wing loading is seen. Also a tighter than 
nominal radius turn is made. 

The lighter vehicle ground track is in Fig. 4.38. 

A higher than nominal turn radius is seen as is a lower 
than nominal landing velocity. 

The high L/D vehicle ground track is shown in Fig. 4.39. 
A larger swing than nominal results from the transition 
ranging guidance responding to the higher ranging 
capability. Also a high landing velocity, 345 fps, shows 
that the flare logic could be improved with increased 
knowledge of L/D. 

The low L/D case, Fig. 4.40 shows almost a direct 
path to the runway with no turn possible. The trajectory 
was nominal up until the flare. There the vehicle fell 
7700 feet short of the runway as the range capability did . 
not exist in the flare. Logic to allow for this possibility 
should be an easy addition. This logic would vary 
BIASl and HF2 as a function of the estimated maximum 
L/D. 
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Fig. 4.3a: Ground Track for 20% Decrease in Weight 
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Fig. 4.39:: Ground Track for 20% Increase in L/D 
























4 . 6 Atmospheric Variations 

Variations were made of i 20% density about the 
nominal density as defined by the 1962 U.S. Standard 
Atmosphere. 

The 20% decrease in density ground track is shown 
in Fig. 4.41. It is very similar to the low weight 
case, Fig. 4.37, as one might expect. Even the landing 
velocities of the two cases compare closely, 306 fps 
versus 299 fps. 

The 20% increased case. Fig. 4.42 is very similar 
to the high mass case. Again, the landing velocities 
of the two cases compare very closely. 
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Fig ,4 . 42 • Ground Track for 20% Increase in Density 












4.7 Navigation Performance 


The effect of navigation errors was determined by a 
simulation combining the steering equations in a closed 
loop manner with an aided inertial navigation system. 

The navigation system is essentially that of ref. [12] 
although reasonably close agreement has been obtained 
with results of ref. [14], 

The system details are in Appendix E. 

The component parts are: 

1. Kearfott KT70 Inertial Navigation Unit 

2. Navigation Filter of ref. [12] 

* with 13 states (2 DME bias estimates) 

• square root formulation 
. non-linear compensation 

m. 

3. Six measurements to update navigation 

1,2 DME 1 & 2 

3. Barometric altimeters 

4. Radar altimeter 

5. ILS localizer 

6. ILS glide slope 

4. A second case using 3 Cubic CR100 type precision 
DME 1 s was considered. 

In the first case, attention was turned to the landing 
site at Cape Kennedy. This landing site and vehicle initial 

conditions are in Table 4-3. 



Table 4-3: Vehicle Initial State and Runway for Cape Kennedy 
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The initial errors are summarized in Table 4—4... 

These are from Ref. ill] as typical of a once-around 
aoort mission using Kearfott equipment. The navigation 
update time was 5 sec. initially, and 1 sec. on final 
approach. A table of significant navigation events is 
given in Table 4-5. 

The ground track for this trajectory is shown in 
Fig. 4.43. The similarity to -the nominal trajectory, 

Fig. 4. 1, is by design. The trajectory parameters are 
shown in Fig. 4.44 through 4.46. 

The altitude error time history is shown in Fig. 4.47. 

Both the actual error and computer one sigma estimate is 
shown. It is seen that there is little improvement 
in this altitude error until the barometric altimeter comes 
on at 215 seconds. A second improvement comes at 940 seconds 
under the action of the radar altimeter. See the expanded 
scale Fig. 4.48. There is but 40 seconds of good altitude 
information before touchdown. But it is enough. 

The East position error time history is shown -in 
Fig. 4.49. The initial value of 40000 meter is reduced 
to about 1000 meters within the first several measurements. 

A second improvement comes with the barometric altimeter event. 
And another improvement comes with the action of the ILS 
localizer, on at 900 sec, Fig. 4.50. Finally, the improvements 
effected by the radar altimeter show their effect in the East 
channel at about 980 sec. 

The North channel errors. Fig. 4,51, are slower to converge 
That is because this is the largely unobservable cross track 
direction. But again, the barometric altimeter helps convergent 
at about 200 seconds. An expanded view of this channel is in 
Fig. 4 . 52. 

The only real effect of these errors is seen in the 
commanded roll angle. Fig. 4.44, where excursions (noise) of 
about 2 degrees is seen. The initial transient as the large 
initial errors are solved is the step from -29 to -15 degrees. 

That the cross track error can be a problem is illustrated 
in Table 4-6. Here are one minute summaries of actual errors 
for the largest allowable initial cross track error (6N=32000) . 
Larger errors cause filter divergence. This is illustrated 
in a failed case. Table 4-7, (SN = -64000). Here, the filter 
settled on a false solution displaced in track because there 
was no information to the contrary. 
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The condition was easily fixed by using the VOR portion 
of the DME for the first 200 seconds. Even though the VOR 
had a bias of 1 degree, this was enough to resolve the 
problem. No amount non-linear filtering will solve this problem. 
More information must be added. An alternate solution would 
be to place the two DI-IE' s further apart, thus, decreasing the 
geometric dilution. 


4.7.1 Precision DME Navigation 

Still another solution would be to use more precise 
DME equipment. Such is the Cubic CRl-00 equipment with a bias 
of 1 meter. 

A ground track of trajectory guided with this equipment 
is shown in Fig. 4.53. The equipment is described in Ref. [12]. 
The three transponders are placed 3 Km to the side 3 Km in front 
and 15 Km in front of the runway. A 30 sec. update is initially 
made with a 5 sec. update on final approach. As can be seen 
from Fig. 4.54 to 4.56, a much lower level of errors results 
than in the preceding case. Were other considerations not 
overriding, such as the cost of equiping alternate runways, 
this equipment would certainly be preferred. 
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Fig. 4.45: Velocity, Altitude and 
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Table 4-^4: Initial Errors 
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Time 


(sec) 

Event 

0 

Start DME measurements (2) 

215 

Start Baro altimeter (H = 30Km) 

905 

Start localizer (H = 2.84m) 

945 

Start radar altimeter (H = 800m) 

960 

Start fast time update (H < 400m) 

975 

Closest approach to DME1 (H = 101.6m) 

981 

Drop DME 2 (elevation less than 1°) 

989 

Drop DMEl (elevation less than 1°) 

1001 

Touch down 


Table 4-5: Time History of Measurement Events 
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Fig. 4. 48 : Altitude Estimation Error & Computed One Sigma Uncertainties 
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East Position Errors (1000m) 



Fig. 4l49:East Position 


Errors & Computed One Sigma Uncertainties 
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East Position Error (m) 



Fig. 4 .-50;East Position Estimation Errors & Computed One Sigma Uncertainties 
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North Position Errors (1000m) 



Fig .-#.51: North Position Errors & 




Position Error (m) 




Velocity error (m/s) 


Position error (m) 

Time 

(sec) 

V E 

V N 

V up 


6 E 


6 up 

0 

.43 

0.0 

0.7 

40000 

-32000 

3000 

60 

23.87 

.11 

1.66 

- 

2655.6 

-22241 

3167.3 

120 

20.75 

- 9.98 

6.20. 

- 

1096.4 

-20858 

4265.3 

180 

16.13 

-38.71 

2.02 

- 

986.91 

-25689 

980.50 

240 

7.13 

18.10 

44.72 

- 

1519.3 

-11268 

831.90 

300 

- 1.04 

45.27 

32.10 

- 

1726.8 

- 6024.6 

1528.6 

360 

-11.28 

28.31 

11.52 

- 

2810.1 

- 4925.0 

819.98 

420 

- 8.30 

11.11 

1.90 

- 

3261.0 

- 3812.5 

108.69 

480 

- 8.63 

3.14 

0.72 

- 

3434.0 

- 3149.4 

- 50.58 

540 

- .05 

3.79 

.83 

- 

2144.2 

- 1338.1 

- 17.09 

600 

5.56 

- 2.78 

.42 

- 

1186.6 

- 672.51 

- 25.12 

660 

4.26 

-10.05 

1.83 

- 

274.68 

j - 725.72 

- 29.47 

720 

- 3.59 

- 7.13 

-1.95 

- 

469.15 

- 900.88 

- 51.75 

780 

- .31 

1.56 

1.26 

- 

500.96 

- 468.48 

- 25.20 

840 

7.10 

.81 

.64 

- 

185.32 

- 104.97 

- 13.76 

900 

2.75‘ 

.21 

.75 

: 

353.80 

- 269.73 

1.24 

960 

2.05 

.98 

.72 

i 

33.98 

8.60 

4.97 

1014 

- 1.25 

- .23 

.04 

- 

9.05 

1.72 

.69 


Table 4-6: Time History of Position and Velocity Errors with 

Large Initial Latitude Error 
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Velocity Error 

(m/s) 

Position error 

(m) 

Time 

(sec) 

V E 

V H 

V 

up 

{ E 


6 up 

0 

0.43 

0.0 

0.7 

'40000 

-64000 

3000 

60 

105.16 

4.99 

1.48 

- 3476 

- 76296 

2809 

120 

106.09 

2.83 

5.85 

2949 

-76522 

3122.8 

180 

96.59 

- 22.40 

2.84 

7418 

-75163 

- 550.27 

240 

93.51 

- 88.94 

- 93.95 

14687 

-88546 

-10745 

300 

19.82 

-202.97 

206.95 

4240.9 

-80166 

10284 

360 

17.58 

-573.22 

19.21 

. 13699 

-132520 

10417 

420 

113.37 

-472.61 

19.77 

20136 

-161300 

4435.9 

480 

371.52 

7.22 . 

57.60 

37796 

-163290 

2186.9 

54 0 

417.00 

199.29 

- 34.83 

78504 

-161940 

1288.7 

600 

773.33 

120.09 

32.85 

136140 

-156420 

534,57 

660 

DJ o . 7 o 

503 ..49 

23.07 

162840 

-134580 

I 

904.24 

720 

426.31 

528.02 

- 6.11 

181350 

-102590 

258.17 

780 

284.50 

693.54 

•36.43 

190770 

-76417 

816.73 

840 

7.92 

626.86 

- 13.03 

186200 

-41012 

587.43 

900 

-189.56 

656.10 

32.77 

177700 

-13577 

•884.22 

960 

-340.35 

521.03 

- 19.25 

159880 

10274 

- 40.25 


Table 4-7: Aided Inertial Navigation with Large Initial Latitude 
Error 
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Time 

(sec) 

Velocity error (m/s) 

V E V N V UP 

Position 

6 E 

error (m) 

5 N 5 up 

0 

.43 

0.0 

0.7 

40000 

-64000.0 

3000 

60 

3.01 

-7 .12 

1.93 

56.45 

- 

2533.1 

2994.3 

120 

1.35 

-*7.63 

-0.93 

47.09 


2245.5 

1889.2 

180 

3.14 

-1.88 

3.72 

443.97 


1922.9 

2884.3 

240 

.34 

1.81 

-1.09 

33.43 

- 

244.72 

301.10 

300 

- .61 

2.15 

-4.03 

10.44 


63.91 

- 122,18 

360 

- .83 

.87 

-3.03 

51.15 


46.94 

- 228.06 

420 

- .94 

- .08 

-1.50 

83.55 

- 

23.85 

- 148.73 

480 * 

- .88 

- .39 

- .68 

97.34 

- 

30.65 

- 134.40 

540 

- .26 

- .20 

- .08 

1.83 


45.08 

- 105.66 

600 

- .15 

- .12 

.17 

11.55 


50.99 

- 73.59 

660 

- .06 

- .005 

.44 

17.28 


37.67 

- 48.31 

720 

.48 

.006 

.67 

19.02 


33.20 

- 30.83 

780 

.60 

- .05 

.64 

47.64 


44.50 

- 24.31 

840 

.70 

.000 

.52 

90.06 


54.45 

- 19.12 

900 

- .33 

- .40 

.50 

45.55 

- 

43.42 

- 13.79 

960 

.48 

.27 

.61 

12.28 

- 

3.60 

15.84 

1007 

* 

- .01 

.38 

.02 

.41 


11.59 

.17 


Table 4-8:Time History of Position and Velocity Errors with 
Large Initial Latitude Error Using VOR/DME. 
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Fig .4; 53: Ground Track for Precision DME 
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Fig. 4,54: ALtitude Error for Precision DME 
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CHAPTER 5 
CONCLUSIONS 


Satisfactory performance has been demonstrated for 
this steering over a large spectrum of conditions. This 
leads to the main conclusion of this report that analytic 
techniques can lead to an effective and flexible Shuttle 
guidance system. The spiral mode may not be most desirable . 

But the spiral approach does put one in the vicinity of 
the airport at maximum altitude, a desirable feature. And 
the circling descent keeps the vehicle in close contact 
with landing strip, also desirable. The near constant 
bank angle mode provides a relatively smoother ride than 
other guidance modes. Wind compensation through center— 
of - curvature control has shown dramatic performance. But 
all these features are away from the main point, which is 
whatever scheme is selected, it should be analytic in 
nature. Each of the flight segments of any of the methods 
studied can be represented by simple formulae. By repre- 
senting the vehicle capabilities with simple formulae and 
then closing guidance loops around these representations, 
a smooth natural type trajectory is almost certain to evolve. 

And new mission constraints will not be precluded as they 
might in some rigid type approach. 

A combined guidance and navigation study has been 
accomplished. These studies showed that satisfactory 
performance is possible with conventional landing aids. 

The convergence of the large initial errors to low 
values is within several time steps and the noise in the 
navigation is seen as small one or two degree roll angle 
excursions . 

But there is a real danger in the use of the navigation 
aids proposed as the standard for this study. Navigation 
divergence can result from large initial errors which push 
the navigation filter into the non-linear and unobservable 
regions. 

In particular, large cross-track errors can cause filter 
divergence. The fixes are obvious, but the trend has appeared. 
Namely, the proposed navigation equipment seems to be tending 
to a bare minimum type system. If this is true, careful attention 
must be given to extreme cases such as the non-linear problem 
that arose in this study. 

A much happier solution to the navigation problems is 
the precision ranging equipment. It is too bad that this 
type system has lost favor because of such things as cost 
and availability at alternate landing sites. 
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CHAPTER 6 


SUGGESTIONS FOR FURTHER WORK 


Some areas of the guidance system were not 
exhaustively studies. The navigation ranging is 
one. The flare logic in the presence of large vehicle 
variations is another. The pilot- guidance system 
interface is another. These as well as all the phases 
can benefit from more work. But in all these cases, 
only minor improvements will possibly result. 

Some larger areas that need work are: 

1. Develop criteria by which to judge the many 
candidate sy s terns . 

2. Use the closed loop guidance system in a 
failure analysis to determine what critical 
parts need be protected and what others to be 
self corrected by the navigation system. For 
example, where would an undetected accelerometer 
failure cause navigation errors that the Kalman 
filter could not correct? 

3. Add the inner loop rotational freedoms. The 
commanded rates are not large but this should 
be checked out. The supersonic and transonic 
lateral stability is always a problem for these 
delta wing gliders. 

4 . Perform man-in-the-loop simulations to further 
determine pilot capabilities and desires. 

It appears that this solution to the Shuttle approach 
and landing problem is fundamentally a good one. It has 
been my experience that the early stages of a program show 
large variations in vehicle parameters, mostly in the adverse 
direction. The latest baseline Shuttle vehicle has considerably 
less flying ability than the 040a. Another large area of work 
would therefore be to: 

1. Adapt the design to these changes or verify that 
these changes are within the guidance capability. 

2. Define performance limits below which the Shuttle 
vehicle may not go. 


Preceding page blank 
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APPENDIX A 

DERIVATION OF TURN ANGLE AS A FUNCTION OF 
INITIAL AND FINAL VELOCITY 


This derivation is essentially the same as that 
of ref. [ 1 ]. An additional term has been added as 
suggested in ref. [ 9 ] to allow for a non negligible 

^ in the equilibrium altitude rate term. Ref. [9 ] 
assumed constant dynamic pressure and the dynamic 
pressure appears explicitly. Here, constant dynamic 
pressure is implied with assumed constant Cj^ & C^ and 
assumed equilibrium. 

If we assume equilibrium flight, the centrifugal 
acceleration is g tan cj> so that the rate of turn is 


d9 _ g tan 4* 
dt V 


(A.l) 


We can differentiate the equation for lift to 
solve for the rate of change of velocity 


dL 

dt 



(V 2 •*£ + 2p V 



(A. 2) 


Since in equilibrium the lift is constant, this 
derivative is equal to zero so we can solve for rate 
of change of velocity 


dV _ V dp 
dt ** 2 p dt 


In an exponential atmosphere, p 
we have 


dp _ 
dt 



and 


dV V£ dH 
dt 2 dt 


(A. 3) 



so 


(A. 4) 


(A-.5) 
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The particle equations are 


• T 

Vy = — cos <j> - g cos y 
m 

In equilibrium, Vy in eq. (A. 7) is negligibly 
so we can write 

0 — — cos <J> - g cos y 


or 


— cos <j> = g 
m 


since cos y = 1 . 

We can solve for tanyusing Eq. (A. 6) & (A. 8) 

dV y D 

j. + 

tan y = - 


1 + sf/£ 

at m 


(p) cos <(> 

We assume tan y = sin y so we can write 


= V sin y = V tan y 

..1 liML 

(^)cos <p 


Using Eq. (A. 8) we can write for D/m 

3 


- S'# -,L 


(fr) cos <f) 


Substitute Eq. (A. 10) & (A. 11) into Eq. (A. 5) 


dV 

dt 


V 2 8 


2 (— )cos <j> 


(1 + 


,L. 

(-)cos <j> 


Solve for 


dV 

dt 


dV 

dt 


V 2 B ' 


2(^)cos <f> 


(1 + 5r> 


SZ) 

dt' 
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Use Eq. (A. 13) to eliminate time from Eq. (A.l) 
<38 _ L sin 4> ^ V 2 Bv 

av-- 2 ’5 -^r <i + Tf> 


- - 2g (fa 2iS_± 
D ev J 


,L. sin 4> 
V 


(A. 14) 


Integrate Eq. (A. 14) 


* 2 - 

^ x UP v t V , 


L V 1 

+ (^)sin <f> In (^p) 


(A. 15) 


Note that constant dynamic pressure is implied so that 

1 „ 2 _ 1 . „ 2 
2 P 1 V 1 “ 2' P 2 V 2 


and 


^1 1 P 2 1 

ln( _l) = Jmi - h 2 > 


(A. 16) 


by using the exponential atmosphere property. Finally we 
can rewrite Eq. (A. 15) by substituting Eq. (A. 16) into 
Eq. (A. 15) 


6 


2 


,L, sin <f> 

B -22 



+ | 8 (j±) sin t (h x - h 2 ) 


(A. 17) 


The first term in Eq. (A. 17) was the form used in 
ref. li]. The second term is the correction due to non 

negligible an( j can account for 15% of the? turn angle as 
at 

described in ref. [9 ]. 
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appendix b 

DETAILED FLOW GRAPHS' 


The succeeding pages show the detailed flow 
graphs for the steering logic. This is for the 
steering program, LNDSTR3 shown in Appendix D.. 

The nomenclature and list of computer constraints 
for this program are in Appendix C. 


- 110 - 


’TERMETRICS INCORPORATED • 701 CONCORD AVENUE • CAMBRIDGE. MASSACHUSETTS 02138 • (617) 661-1840 




-?T7\e.T furPAUfce 































US' c 



VUI2»jl , jc.Cu?0- ' 


tffti. = j,-_( 121) r,ct> 7 ? 

Kttirr^ /«.('. v) - > c ^ 

byrr 4>tow- -^'0 

vltf/- /*.('}<*) - 7 


\A/!> : 6* (96) - Sooc? 


Ca-TT- 6 'rt.C 9 ' ) - .'= 


VF ’ (*('•', \* 1°™ 

Ai 7 !.’- iv.'V'.)-. l-J, 
PQ>'/!- >-T 

COj ‘ f * c ( r : J s 'tr 











- 114 - 





















HF loooo 

- A5~ 




^esTrCi, OPbTZDU 






tit a > . u<co}/$ 

'L/^fT^a DCTW/JK) O^D/z. 


X>gc 30 * jkc.(j 5 D)= .?< 


0): u>i • VKTl 9 *=■ 

P 3 ** (VaV^* uof i 

_ JL 

H Sr^(v*i/errt) 

Oi = DD 4 P'S 
tflATr V'Oeri. 


< TVeM LC QIC,. 





J>dcl* -wtl* Stff***) 
w:u .- peuTc^r+pJ^u- 
rcrrr i 


ECU-= ATAtJ /jgk ^ A ei3> -p£^r 

»T i 

LipnrRTuL it? £Lr 


* , 4 *- (H ?; ) «= 3 $ 


1 

K 

^ tm^r- ff 

>^H 

ft«u.« k«> no + kw vu\r 



Mi 

cn 






!/*^n4'' - 


1 





Ko\)’ 

JU>t.Oi57 •= 

— ♦ oi7 


</ /K< Vl<y v 

1 j 

1 L6f 

-il 

f VLfr - 

Let. (110 T 

9>6 


M 

Lt<C •< Ptc/,o^iA ,1 ^'^If 
Lim/T tv IpO w 
WT <4 )• . -n*A AJ >0 



Htt= nri+ («m-VKPn. l )/?^rt 


LK> r JUc (l 4 T)= *> 



Ha? Ut(\yi>)-si> 

<Yfls Cc<(lM) « i O 
KPn- £rt(in) "tS 


fc>i.(!lt)r < 2 w» 


v- 

$KL~ Wi e LAW 



■* 3 i* ixCME') - 70 
VCO*- A»:.Cu^)r 4i5* 
tf-tiH- £fc(di)^ -cc 3 
“ ietimy .r 


4tRL rCstf'* 5 ?) * 6 
£C1>\ - Xee(H5) = 1.4- 
KA'_F = £*'c?-'») *=■ 2 
JCJ’-vi'- jCre.O rt ^i *■02 

4-.- ”t.\' w*. v !+•)» 4 
ALFm -"-£*• U4-)- I 1 ? 


FU.MIC IA*JJ 



LnJp 7 )r 27 CO 

RDmiJ’Art^ r -zT 

POTMX ” jfpc-0’^) " -tS" 

RCT?’ - i<«,(llC>) r -4- 

<*?*,- Uz. (<P!>) =6 

WTTz." Ui. ( ,0 f) c 7.0 


= /cc.(i?lj- to 






















- 117 - 


ourror* errueN 



CCL? 


Kmr- A* 


TT? 

tfC trif.T\)& 














Appendix C: Nomenclature for Computer Listings 

VARIABLES FOR STEERING PROGRAM 


NAME 

UNITS 

AMAG 

FPS 

alph 

DEG 

ALFO'.D 

DEG 

A7I-M- 

DEG 

CA/R 

N.D. 

SAZ» 

N.D. 

CTH 

N.D. 

STH 

N.D. 

CRL1 

N.D. 

CTEST 

N.D. 

dale 

DEG 

DA7 

DEG 

0 

FT 

D EL 

FT 

ORL 

DEG 

DEI 1 

RAD 

DEL2 

FT 

DEL3 

RAD 

DH1 

FT 

DRL 

DEG 

DT 

SEC 

DT T c 

DEG 

DTH 

OEG 

DTH1 

DEG 

DTHO 

DEG 

D2 

FT 

03 

FT 

FCTRl 

ND 

GS 

FPS - 

H 

FT 

H : 1 

FT 

HT 

FT 

HS 

FT 

HI 

FT 

IND 

N.D. 

IN02 

N.D. 

IND3 

N.D. 

IND 4 

N.D. 

INDIO 

N.D. 

TNDT 

N.D. 

INDT1 

N.D. 

r- 

N.D. 

DFL 

RAD 

IH 

N.D. 

LG 

N.D. 

LAT 

DEG 


description 


TOTAL ACCELERATION 
COM-.ANDED ANGLE OF AT ACK 
OLD COM* ANDEO ANGLE OF AT ACK 

azimuth 

COS OF RUNWAY AZIMUTH 

SIN OF RUNWAY AZIMUTH 
COS TURN ANGLE 
SIN OF JIJRN ANGLE 
COS OF RL2 

INDICATOR OF CENTER CONTROL 

CHANGE in alpha 
A7IMIJTH ER ‘0 ' 

LATERAL DISPLACEMENT 
CALCULATED LEAD ANGLE IN S**TURN 
ROL ANGLE INCREMENT 

TURN ER % OR FOR 180 DEG CASE 

DISTANCE BETWE. N DESIRED AND TRUE CENTER 

LEAD ANGLE IN s-turn 

CALCULATED ALTITUDE LOS , IN TURN 

ROL - ANGLE INCREMENT 

TIME INCREMENT 

TURN ANGLE (DTHO-DTHl) 

PREDICTED TURN ANGLE 
DESIRED TURN ANGLE 
PREVIOUS VALUE OF DTH1 

LATERAL DISPLACEMENT FROM RUNWAY 
LATERAL DISP WITH NOMINAL ROL. 

FACTOR IN TURN ANGLE CALC 
CONVERSION FROM 6»S TO FPS'- 
ALTITUDE 

ALTITUDE IN TURN ANGLE CALC 
ALTITUDE TO START FIRST FLARE 
SCALE HEIGHT 

ALTITUDE IN 180 DEG TURN CALC 

INDICATOR FOR MODE 

INDICATOR TO RETARGET FINAL PHASE 

INDICATOR TO START S TURN 
INDICATOR FOR S TURN 
INDICATOR FOR FLARE 

INDICATOR FOR ROLL REVERSAL IN TRANSITION 
INDICATOR FOR RETARGET AFTER TRANSITION 

SELECTED MODE 
AZIMUTH ER'-'OR 
INDEX IN DTH1 CALC 
LANDING GEAR INDICATOR 
TARGET LATlvtjDE 


I.NGT 
LOC ' 
OM 


DEG TARGET LONGITUDE 

N.D. COM '.ANDED L/D 

RAD/SEC TURN RATE 


1.1 3 - 





RTOGO 

FT 

RANGE-TO-GO 

RTOG'Ul 

FT 

NQMI JAL RANGE TO GO 

R nT^EF 

FPS 

REF R"NC>' ALT! 1 IDE RATE 

ROT 

FPS 

ALT I' i IDE RATE 

PTQ 

DEG 

CONVERT 'ADIaNS TO OEGRE S 


FT 

FARTH RADIUS 

P')»_ 

DEG 

CO 'I ANDE • ROL 

R 

FT 

FINAL TURN i ADI is 

R 

FT 

IN'T AL TORN AO Is iS 

RL7 

DEG 

COM A NOE > ROL IN TRANSIT* ON 

RL2 

DEG 

COM' ANOF’< ROL CONSIGNED) 

P ^ 

FT 

DISTANCE DETWr: N DESIRED AND ACTUAL 

S* 

FT/RAD 

CONVERT VAO TO Pf 

SB3 

% 

SPE D DRAKE T . <M 

SR 

% 

COM ANDE ; SP.E O DRAKE 

TAMRLl 

N.D. 

TA >J ( RL1 ) 

TANR1 

* NJ. D. 

T AN ( RL<? ) 

TANR1.5 

N,D* 

TAN(PL-i) 

TFS r 

n.d. 

D-STRED SIGN OF R‘)l. 

TFS r l 

N. D. 

DESIRED SIGN OF ROL 

\/ A 

c-pc; 

*/ c l.OC ITY MAGNITUDE 

VAX 

FPS 

-^I.AtI'/E > LOCITV MAGNITUDE 

x/1 

FPS 

VELOCITY A" URN STA i ’ 

VAl 

FPS 

f-'OUTL'MR* iM '/Fi.oCTTv a t ’Li? 

'/IX 

PPS 

TNE°T JAl. VELOCITY MAGNITUDE 

''1 

PPS 

VI IN r CALC < 

vlat 

crpS 

LATERAL VELOCITY 

'J? 

FPS 

/EL OCTT Y A:- ENO O' 7 TURN 

V2 

FPS 

>/2 IN R CAL ' 

VRSR 

FDC^ 

velocity . N so calc 

XI 

- 

DUM Y VARIABLE 

x2 

— 

OUM Y '/ARI-vQLE 

VECTOR 

variables 


R 

ft 

POSITION VECTOR 

V 

FPS 

VELOCITY VECTOR 

URT2 

N*f). 

PROJECTED TARGET 'JNTT VECTOR 

UC1 

N.D. 

UNIT VECTOR AT DESIRED CENTER 

UNI1 

N. D * 

UNIT NORMAL TO PLANE OF VEHICLE AND 

iJR 

N.Dt 

UNIT POSIT ON VECTOR 

UC2 

N. D , 

UNIT NORMAL TO PLANE OF ijC AND UC1 

OE 

N.D. 

UNIT VECTOR EAST 

DEO 

N. D. 

JN T T'AL UNIT VECTOR EAST 

UET 

N* D . 

UNIT VECTOR EAST A ; *WGET 

IJNI3 

N.D. 

UR X UN 11 (CROS PRODUCT) 

*JRT3 

N.D. 

IN* T ' AL HIASEO UNIT ARGET VECTOR 

HL 

N.D. 

UNIT NORMAL TO PLANE OF ■> AND V 

iJhiT 

N.D. 

UNIT VECTOR NORTH A) .ARGET 

URT1 

N.D. 

UF X IJRT (CROS PRODUCT) 

UF 

N.D. 

UNIT VECTOR ALONG P IELD 

UC 

N.D. 

UNIT VECTOR AT INSTA 1 *' ’EONS CENTER 

IJZ 

N.D, 

UNIT VECTOR NORTH 


T/WGfcV 


END OF DAT A 
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49 

51) 

7 b 

r* 

Mil 

HI 

n2 

•H 

v.» 

b*> 

8 7 
8 

b tJ 

wi! 

93 

92 

93 

94 

9b 

9b 

97 

9ti 

9 
1U 

12 ’ 

102 

103 

104 

105 

10b 
107 
1U8 
109 
1 0 

15 
1 '2 
113 

1 1 4 
l’b 

1 ! 6 
157 
118 
1 i9 
12U 

121 

12 > 
• 123 

1 24 
123 


INDEX FOR CONTROL VARIABLES FOR LNUST3 


1/31/7-5 


VAR 

ALFC 
nl AS 
1LT 
.<Df *F 

•\IJMPT 

1N022 

kbiast 

KHI&S.J 

tUEL 

.191 

I ■1015 
•)A«.r -1 

MCI 1 
% J AST 

-) -4. 

9L8 

, wL7Mf| 

ACTR 

'/F 

DELTR 

VLM3 

BIASFL 

RDTMX 

ROTMN 

KSB1 

HSB 

RL2MN 

KCN 

KBLND 

HS1 

FCTR2 
KCL 
HT 2 
OF 3 
RQT3 

KSB 

KSBH 

VRO 

KV-' 

KO • 

VLG 
I LG 
SB 
ISB1 
ALF7 

KTR 

VLM2 

VLM1 

kturn 

IN012 


UNITS 

DEG 

N.M. 

N.D. 

FPS 

N.D. 

N.D. 

N.D. 

N.D. 

N.D. 

N.D. 

N.D. 

DEG 

DEG 

N.M. 

NAD/S 

DEG 
DEG 
FPS . 
FPS 
RAD 

FPS 

FT 

FPS 

FPS 

N.D. 

FT 

DEG 

N.D. 

N.D. 

FT 

N.D. 

N.D. 

FT 

FPS - 
FPS 

N.D* 

1/SEC 

FPS 

N.D. 

N.D. 


VALUE 

6 

10 

1. 

-250 

.0 2 

1. 

0 . 

0. 

0. 

n. 

1. , 

.3 

4b 
30 
.01 ■ 

30 

lb 

lb 

10 

• b 

30 i 
120 i 
-15 
-25 
0 

1200 

15 

• 0003 
0 

310 U 

1. 

1. 

70 

6 

-4 

• t> 

.003 

425 

-.4 

-.017 


FPS 

N.D. 

PERCENT 

N.D. 

DEG 

DEG/FPS 

FPS 

FPS 

N.D. 

N.D. 


5t>0 

1 

20 

1 

7 

.0075 
40U:» 
120 • 

1 

1 


DESCRIPTION 

NOM ALPH - FINAL PHASE 

BIAS ON TARGET - BEFORE FINAL PHASE 

INDICATOR FOR LEFT 1URN LOGIC 

REF ROOT DURING TRANSITION, 

DAMPING « AIN DURING TRANSITION 

IND. FOR COMPENSATING ROLL TRANSIENT •! ! »*H \L°-< 
GAIN TO BIAS FINAL PHASE ALPHA FOR «IND 
GAIN TO BIAS TARGET COMPENSATI G PUR -l '-'US 
INDICATOR FOR USING RELATIVE '.LOCliY 
TURN COUNTER 

INDICATOR FOR CENTER OF CURVATURE CO'-HML 
LIM-.T ALPHa RATE (SIMULATING vJTO^JL » T > 

LIMIT ROL' IN TRANSITION 

BIAS ON TARGET BEFORE FINAL ‘-’H^St: 

MAX AZIMUTH RATE SIJPERSONICAL * 


MAX ROLL IN TRANSIT {ON 

MIN ROL- IN TRANSITION 

NOM AC-,’ EL DURING TRANSIT 'ON 

FINAL VEL FOR TRANSITION RA -GE P^c'HlT. ON 

limit HEADING ER 0 j IN TRANSIT 

VELOCITY TO START AZIMUTH PATP Ll^ i* 

BIAS ON TARGET DURING FIRST FLARE 
MAX ALT RATE IN FLARE 
MIN ALT RATE IN FLARE 
NOT USED 

MAX ALT FOR SPEED BRAKE 
MIN VALUE OF RL2 
GAIN FOR CURVATURE CONTROL 
BLENDING GAIN FOR ALPHA 
SCALE HEIGHT IN TURN CALC 

FACTOR IN RR CALC. 

FACTOR IN DTH CALC. 

ALT THRESH FOR 2ND FLARE 
INCREMENTAL G FOR 2ND t-LARE 
ALT RATE AFTER 2ND FLARE 

SPEFD BRAKE GAIN 

FACTOR IN NOM VEL FOR S.B. 

NOM VEL IN S.B. CALC 

GAIN IN FINAL PHASE LAT CONTROL 

GAIN IN FINAL PHASE LAT CONTROL 

VEL TO DEPLOY LANDING GEAR 
IN FOR LANDING GEAR 
SPEFD BRAKE SETTING 
IND FOR SPEED BRAKE 
ALPHA AT END OF PITCH OVER 

PITCH OVER RATE 

VELOCITY TO START PITCH OVER 

VELOCITY TO END PITCH OVER 

GAIN ON DTH-DTH1 -120- 

IND OF TANR11 IN R»-. CALC 



12b 

ROL 4 

DEG 

20 

12 7 

1* ISM 

N.D. 

1 

120 

IND8 

N.D. 

0 

129 

L'*b 

N.D. 

-1. 

IbU 

Hi AS1 

N. M. 

1.0 

lb* 

RP2 

FPSS 

10 

lb2 

^OT2 

FPS 

— 1 0 


H}-2 

FT 

bO 

lb4 

Al.r 

DEG 

0 

1 b-> 

J TM--J 

DEG 

4b 

1 bo 

UTHLMl 

DEG 

90 

1 b7 

INU6 

N.D. 

1 . 

).bO 

IND7 

N.D. 

1. 

Lb9 


— 

— 

1 a.j 
1A 

K0MP1 

alf^n 

DEG 

.02 

4 

142 

alfmx 

DEG 

lb 

14b 

R TGLM 

FT 

60 ‘ 

t A 

!)THC 

DEG 

20 • 

1*+D 

UL • A 

1/DEG 

1.4 

140 

KALF 

N.D. 

2 

147 

LDO 

N.D. 

b 

148 

RLb 

DEG 

bO 

14V 

DTBIAS 

DEG 

20 

IbO 

UELb 

RAD 

.7b 

lbi 

D2LIM 

FT 

20Uu 

102 

DEL-IM 

FT 

200 i 

lbb 

LM2 

N.D. 

2 

lb4 

OTHLM 

DEG 

30 

lb 

K.OMP 

N.D. 

.01 

lbb 

ALPH 

DEG 

29 

lb7 

TMAG 

— 

— 

lbO 

MOOT 

- 

— 

lb9 

M 

SLUGS 

43b0 

lfoU 

LOD 

N.D. 

7 

lbl 

AZR 

DEG 

75 

162 

RL1 

DEG 

30 

163 

- 

- 

— 

164 

FTR 

N.D. 

.9 

16b 

HF 

FT 

lOOu 




LIMIT ROL TO CHANGE 4LPH I 'I El AL PHASE 
NOT USED 

INDICATOR FOR USING ROL 4 
LIMiT ER O , ON C COOT* )L 

bias on targe r jn F) «al ph-»se 

INCREMENTAL G FOR FLAhE 
NOMINAL SINK RATE AFTER *-LARE 
FLARE ALTITUDE AFTER LA i‘~. 

NOM ALPHA IF - > A'.T 

BIAS TURN ANGLE -fO b‘AR r FINAL PHASE 

BIAS TO STAR i STRAIGHT SEGMENT 

NOT USED INDICATOR FO < U1 'ECT 0 r-INAL 

INDICATOR for BIAS'NG ALPHA-FINAL PHASE 

DAMPING >'-A JN-F J -’aL pH \SE 

MIN ALPHA 
MAX ALPHA 

RANGE TO END OF SET 'ARGE f CALC 
MIN TUR J A-'iOLE FOR C ,ONTR-»L 
CHANGE OF L/D .VI TH ALPHA 

GAIN O'! DELTA L/U-FINAL PHASE 
NORMINAL L/D-FINAL PHASE 
NOMINAL L/D-FINAL PHASE 
BIAS ON PREDICTED TURN A jGLc 
NOMINAL ROL. IN b-TURN LOGIC 

MIN ER: :0 ■ FOR S“TURN 
MIN ERF OR FOR CF CONTROL 
LIMIT ER -'OR FOR CC CONTROL 
BIAS TO END STRAIGHT SEGMENT 
DAMPING GAIN-PHASE b 

NOMINAL ALPHA-BEFORE FINAL PHASE 
NOT USED 
NOT USED 
VEHICLE MASS 

NOMINAL L/D-BEFORE FINAL PHASE 

RUNWAY AZIMUTH 
NOMINAL ROLL 

FACTOR IN DTH CALCULATION 
ALTITUDE TO START FINAL PHASE 


END OF DATA 



APPENDIX D 


COMPUTER LISTINGS 


The Fortran computer listings by which these 
trajectories were generated are included here. The 
main program is EQM0T1 . It calls in turn: 


INIT 

A 

ATMOS 

LNDSTR 

AERO 


The initialization program 

The input data file 

The 1962 U.S. Standard Atmosphere 

The steering program 

The aerodynamic subroutine 
AERDY4 contains the 040a vehicle 


Sub. Fort contains the following subroutines: 


DIFEQ 

ATMOS 

AXB 

UNITV 

VSUB, VAD 

VECSCL 


4th Order Range Kutta Gill integratter 
The 1962 U.S. Standard Atmosphere 
Vector cross product routine 
Routine to unitizing a yector 
Vector subtraction and addition 
Vector times scalar routine 


Subl. Fort contains these subroutines: 


SIND, COSD 
LALONG 

DANG 


sine cosine of angles in deg. 

converts unit vector to latitude & 
longitude 

calculates angle between two vectors 
by cross product 
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1 

(-■ 

K> 

U) 

I 


2.1069C56F 07 0.0 0.0 

-0.0 — - — 0-.0 7.9 71CA//t-01 ■ 

1 .OCCOOCOE CC 0.0 0.0 

2. 5972724E-05 2.0000000F 00 0.0 

t> . CCCCOGOF CO l.COOOOCOF 06 0.0 „ • 

1 . 5 OOOOOOF 00 5 .OOOOOOOE 00 3.0C0000Ot 01 
0.0 1. OOOOOOOE 00 1.2C0C0U0E 03 

-0;0 — l . 5000000E ^0 o.O CCOOOOt 00- 

4 . CCCOOOOE OC 2.2 7V9908E 01 2 .5 CrOOuut-Ul 
1 . 2 C 000 COE 04 0.0 7.5CC00UuE Oi 

-5.0CCOOOOE 01 2.00000007-02 0.0 
2 . OOOOOCOE 02 2.C0000CC C 03 ! .OCCOOOuc. Ou 
2. 3090000F 04 2.250000CE 04 3.5C00O00E 0-t 
• 3 .469999C c -02>-5.5509996E-02-9.0'*79n : ,vr -02- 
-2 .37 19 04 5E-C 1-3 .3049 9 95E-0 1-3.6 04 99-lot -ul - 
•2.02C999CF 00-3 . 353008?= 00 l.CCTCOOOt: 00- 
l.OCOCOCCC OG 0.0 0.0 

1 . OOOOOOOE 00 5 . OOQOOOOF-O 1 4.5CCO000l 01 
3 . OOOCOOOE 01 1 . 5000000E 01 1.3COOOOOE 01 
-3.0000000F 03 "1.20000Q0E 03-1 . 5 OCOOOOl 01- 
1.2000C00E G4 1 . ‘300 OPC<'>' : 01 2 . 9 99 9*-'9ofc -v ♦ 

1 . OCO r, OCOF OC 1 • 0000000 r 00 7.0 CCOOOut 01 
b . OOCOOCOE— 01 2.9y<,955flr-03 4.2 5CC0U0L. U2- 
b . 5 n COCOOE 02 1. OOOOOOOE 00 2-OCCOOOOt 01 
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2 . OCOOOCCF CO 5 .OOOOOOOE 00 3.0rC0O00e oi 
2 . OOCOOCOF 03 0.0 2.0C90000E '00 

-2 .90000CCF- 01-0.0 -*-0.0 

7 . 5COGGOOF Cl 3.0000000*= 01 0.0 


2. OOOOOOOE 

02 

6. 0:> i<i 7^0£ 03 

00009005 

/. 0710 6 7/C- 

-or 

L.QOOOUUUl 00 

cc :c ) c 1 3 

9 . v e2*rb t)4 c- 

-ol 

v j ioOz-02 

00.) ) ) 115 

1 . uuO JOuut 

0 0 

<+*:>OOOGJOc 01 

000 >>020 

2.000jOuUL 

0 1 

t>. OUG J U UOc 01 

0 00 00025 

1.200000UE 

03 

l.OOOOOOJc 00 

00900030 

1 . oOGOOuuE 

00 

1. OULU 000:: JO 

00000335 

2 . o^yWoGw 

00 

i.OOOOOOOfc 00 

QOU-JOOtO 

D. OOOOOOOE 

00 

2-000J00Jr 00 

QUO 700 O 

t>. OOOUUOOL 

Ou 

l.OOCOOOOc 01 

000 300-j3 

0.0 


1 . 000 J UOuE O’J 

09090095 

1. OOOOOOOE 

00 

2*0^900UE 04 

009 )10u0 

0.0 


0*0 

00090065 

1 .4 099902 E 1 

-01- 


90 .’390/0 

4 . 9b 59 9 ) 7 l 

-01- 

■ 0 . v Soz.c-0 i 

00 )0 10 / j 

2 . o 0 0-J(- \ 

„02 

2 « iju J .*) ,) J W T. J L- 

09 J»V)3 ", ') 

1 . 0000 OOG c 

uO 

1*OoJ*3uO>jli J 0 

000 J JlOu 

3.O0000O0C 

01 

1 * 000 JO JO? 00 

00000090 

1 .OOOOuOuF 

03 

5* OuOoOOoc -0 1 

000 10095 

2.3 O0000 OF. 

01 

0.0 

00 09 1100 
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3 . LO 0 J J 0 ot: O'-t 

J J « J 3 0 1 0 5 

o . 0 00 J 0 OU C 

00- 

•^•oOJuUOOc ^ J 

00 Ml 113 

J.')vv'/v JiiL 1 

-01- 

■ 1 • wL ) j V ^ J ii ^ J 2 

09 On 0 l U> 

1.0 000 0002 

00 

7.oooouuor oo 

00 POO 120 

l.UOUuOO'Jc 

00 

l.OOOJOOOz 00 

00 10012 5 

1. 0000000-' 

Ou 

I.yy*;wd3u ou 

IP 1 39 

u . jO J J *J Jvi. 

00 

*+ * b u \ t ) U J u — *J L 

) ' i ; i ) > j 

3.1 bi2 9v<~>fc 

03 

l . 0 0 J J vo «r — i J<! 

>.,■'•>101 * j 

2.0 OOJ 0 OU L 

02 

1# J /l! 0 0 

OG o 10 1 c 5 

2. JOOuOuOn 

01 

7 . bOO JOuuE-U L 

10 iJ'JLsO 

3 .OOOOOOOE 

03 

v - 9 v v v 9 '♦ -0 3 

090 1015 5 

4 • 3 bOuOOu E 

0 3 

7 . 00 JO 0 'JO c 0 0 

009 391 ->9 

d. 994999^,^ 

-01 

1 • UUoO u OOz J** 

000 )Ui>5 


File A: Input Data for Nominal Case 
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LNt/iTFi 


DATE '= 730 78 


15/55/44 


CC02 
' CC03 
0C04 
0005 
“0006" 
0007 
0003 
"0009 ’ 


NJ 


0010 


* 00 11 “ 
0012 
0013 
*00l<» 
0015 
OC 1 6 

0017 

0018 
0019 
0C2C- 
0C>! 
0022 
0023- 

0024 

0025 
- 0026 “ 

CC27 


SCR" CO 1 IM 1 lni>5Trt(K,V,UcO, VAX.VlX.AtlM.rWRULL, ALPH.A^AG, 

VI, V2 ,01t- ,rTHl,l>33 ,1N1 ,LATI, Lu!m 37, Aif-'l.^ruGo.^r-l, ALF1,0T,LG,'SU> 
Implicit 3 = A( U-m) 

INTFGFC I , Th, 11, ll,u» lilAKI 
CAT4 ISTA = T ,9 Rl/0,0./ 

CTMcrsirf! o(3),V(3) ,A( los) , JtO( 3 I ,tJL (3) , A/(3 ) , TL,-l{ 3) ,UNI U 3) 

C.J «fA SI~m USt (3) ,OtlAl ,ukT (3) ,IJ!;T (3) , JF ( 3 1 ,0* 12 ( 3) , 'J* ( 3) ,UNI3(3T 
cr-'cwsrrv L°Ti(5i,uCi3) ,oci ( 3 ),ut 2(3), u*rj(.j> ,b( 12 ) 

0JW4SKN L51 ( 3) , IE K 1(31 , Ufc ( 3 1 

EOtJT VAL ?N‘C9 (A<l4l ,Okri, (AC/1 ,U21 , (A(201 ,KL31 , (A<2Q),PL4), 

( A (83 1 .KB IAS*)}, ( A(o2),K.olASr 1 , (AC 781 ,1U1 , 

(A (80) ,K0*'FT) f (AC 7V1 , RO TRET 1 , ( A( o l ) , IND22 1 , 

<4(95),n=LTP),<A£9t),*F> , ( A 1 93 ) , ACT* ) , ( A ( 92 ) » AL 7.lN) , 

(A(37),9M.FL v >,(A(oo) ? 1(1013) , (»(o5) , In 1 ) , ( A( cN ) , I PEL > , 

(A( 91),t t C ) , ( A ( VO ) , UMcil 1 , ( Ato9> ,o!Ai> 1) , ( 4 (tie) ,KLTl 1 , 
(A(99)» 0 C'T*'\ , )»(Al9u) »kOTKa 1 , l A< 97 ) »i)IA:»FL)»(A(9{>) »VLM3 ) » ' " ‘ 

I A( 103) ,KCM,(A( 102) »KL2MN 1 , (rt ( 101 ),Hbu) , ( A (100 ) , KS3 1 1 , 

(A(10 7),KCL), (A(lOo) ,FCTR2), ( A( 103), USD , (A( 104) .KtJLND), 

(All 1 1 1 ,KS a ), ( A( 110) ,R0T3) , (A( 109) ,GF3) , (AI103) ,i(TT2) , 

(A( 115) ,Kf>C)»lAtlU) ,KVV)r ( A ( IIj ) , Vh 0 ) , ( A( 112) ,A*dH ) , 

( A(1 14) , I 5P1), U( lia) ,ii>l) , (A (117) ,1 Lo) , ( A(llo) ,VLG) . 

(A( 122) ,V1.vu ,(A( U2), VLM2) , ( A ( 12 1 ) , KTR ) , ( A ( 120) , ALF7) , 

(A( 126) ,p.p| 1.4) , t At lu), LUiNiO ) , ( A ( 123} , (ND1 2 ) , ( A( 12t ) ,KTURM) , 

(A( 130 ,R IASI), (A{ 124) ,L.>13) , ( A(128), IN03 > , (A ( 12 7) , IN09), 

* (All 34),. M FF1, (A( 133),HF2) , ( A (132) ,RJT2) , ( A( 1 3 1 ) , CF 2 ) , 

* (A( 137) , I ACM, (A(1J8) , ING7 ) , ( A(13ol , JThL M .l), (AC 135) , DT'IN), 

* (A( l61l,AZ R ),(Atl<,2)*i\Ll)i (A(loO) , LOO) , ( A ( 164 ) * r f A ) , ( A ( 165 ) , HF) , 
*“ (A(19),PtJTP), ( A ( 21 ) , PCTR ) » ( A ( 4 9 ) ,ALFC) , ( A (50 J , t> I AS ) 

tOl/'VALEKCF (A(l5t>) ,ALFC1) , 

* <A(154) ,PTH«) f l A( 1»») »KUMP) , < A< 15 3) , LM2 ) , < A ( 1 52 ) , OEIL 1'1 ) , 

* t A( 150) .DEL 30) , (A( 131) ,021 I M) , ( A ( 149 ) , OTb I AS ) , ( A ( 143 ) , RL5) ," 

* ( A( 144) .OTKC.) , Ull<t3) ,OLL)Jm) , ( A { 14o) ,KALF ) , ( A( 147) , LOO) , 

* (A( 143) .am*') ,i A( A 42) , ALFMX) , (At 1*,1 ) » ALFrlN ) , ( A ( 140 ) , KDMP1 ) 

DATA LG/1.,0.,0./ ’ 

IF ( T $T A° T . EQ.l ) GU TO 111 
0TH=36I . 

UT^ 1 1 = 26?.* 

CTHr>=rT)-l 
ISTAPT=1 
— I MO T 3 0. 

P 1 3 3 . 14159265 
PTD=ieo./PT 
GS = 32 .2 

A-'A'UOS 
FS=25C00. 

.... C ^S’LC= 1 . 

SF=»TC»60.*60tjO. 

PF=205C290C./60b0. 

- * - 

FCT° 1 3 1 . 



FGPTPAN IV G LEVEL 21 


LNOSIK 


DATE 


OC?8 — FFAtV (7'li2U'A ' ’ " 

0C2« CALL AXBCL7,l)tvl ,TE,-i) 

0030 CALL L-HtTV(TFM,oer,)U f XlJ 

0C31 CALL AXO(L°T*Ubl ,UNI ) 

0032 TAN31 1=S]NC (Sf*) /CO^OC RL4) 

0C33 TA MO 1 l=T AN^Ll 

0034 T7N->t.5 = SI'‘0( c L5J/CG5D<RL5) 

0035 $A7«* = SI!*n<A2 :, l 

co it CPU=rrsc(PLAi 

0037 C A Z° = CCSO ( A Z° ) • 

0038 ASSIGN 1 TC 11 

0039 IF(H.Lr.HC) ASSIGN TO II 

0040 1 C (K.LT.HF) INI>=* 

OC-,1 0 n 11 1 = 1,3 

0042 ’ 11 UF( M= (UNT ( T)*CAZR+UbT( I )*SA1R) 

0043 CALL AX fl( UF i UKT 1 ) 

0044 DC 1116 1=1*3 

0045 1116 TEN'( I )=U°T 1 1 1-JF ( I } *oI AS/RE 

0046- — C.'l L LN'ITV(TFH,lMf3./U,xn- - 

CC4 7 21 F(1 D 3A I (5T14.TI 

0048 0^ 1110 1=1,3 

0049 1119 U' 5 T2{ I )=L , ’’T3( 1 ) ■ 

0050 IFIVI X.LT.VLM11 GU 113 111 

0051 DO 1222 1=1,3 

■00 52 1-22-2 TE'M I )=U= 7 Ul-OM 1)*BI AST/RE 

0053 CALL LN'ITVITFH.GKrx.Xl.Xl) 

OC 54 111 CALL AX3I v,4, Ic.’ll 

■ 0065 CALL- UNIT VCT6HniL»M»Xll — 

0056 CALL LN!TV(°,Un,Xi,Xl) 

0C57 CALL AX8(L'2,lW,rbMJ 

0058 CALL LfmV(TFH,Ot,Al,Xl) 

0055 CALL ACCT!MUF,0L ,X2J 

OOoC V A= V I X 

0061 IFHPEL.SC.l.) VA=VAX -- 

0062 VA 1 = V A 

0063 AZI M =A3CCSIX?)«AI J 

.0064 CA| L AOrtf»( V,o:,/i) • • - — • 

0065 IMXl.LT.C.) AAl.".=sbO.-AZIM 

CC66 IP(!NC,FC.4.ANi«.lAD2.NE.l.l GU 70 4*t4o 

0067 -4448 CALL AXP t L' c T? *n , I ci-' ( 

0068 CALL LN‘ IT V ( TFM,u.Nl 1 » Xl , XI) 

0C65 CALL AXBIUP.UnH.UNIJI 

. 00 70— CALL . ACGTHIL'N’l J,UL* ARGJ 

CC71 | TFST *SI'7N(l, l A»vGi 

0072 ‘ CALL «r)'T , >(V.JK,K!JT ) 

0073 CALL AX :M U c T? ,G!' , ThM) - - - 

0074 .... CALL LNITV(TEM,UNU,Al,Xll 

0075 CM L /‘fTTSIllMli.ULfXn 

-0076 IFlXl„r,T. 1. ) Xi = U--~ 

CC/7 • CEL = A 3 C C S (XI ) 


73078 


15/56/f4 


PAGE 0002 
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LND6TR 


DATE * 73078* 


15/55/44 


PAGE 0003 


pr 78 
00/9 
CC8C 
ooat 

QC*? 

co a 3 

— n-;v 
0*. 5 
0CT6 
0027 “ 
0088 * 
0CJ9 

— Ci rc r “ 

CO u 
0CV2 
0003 
CC94 
C095 

— 0006 — 
COS 7‘ 
0^98 

* 0C99 *“ 
OICO 
0101 

— 0102 

01^3 

0 1 C A 
. C 1 0 5 

H 0106 

N> 0107 
Oi :g — 
1 n:r> 
QUO 

out 

0112 

0113 

- - nil'*™ ■ 

01 15 

0114 

'*■ 0117 — 

oua 

0119 

01 2* — 

o:n 

G17? 

... 0l?3 

0124 

0175 

..... Q12t . 
0127 


cmV^-CA^ctup, 6 ri2)*sf 

GO TP n.<U2.3,*> 
inVIX.CT .VLM1) 00 TU 5558 

ifi n’rrurc.u) go to 555 * 

DO 5553 T ” U 3 

5653 T£*M 1 )=lPT3U)-(N8lAir^tVIX-VAX)/RE}«UF( I ) 

~ CM l L\T T V ( T r n t Ur.T2 r Xl > XI) ” 

CM L LMTV(iJM^,uKMfXl,Xl) 

r >p *u=i. _ _ 

GP TP 4446 " * 

5554 RPU = &L 1 * T£ St 

nv = CS*$UC<M1 J/ICUaOlRLl )*VA> 

T*=unv*nn .CT.Db L ) RJLL=AT AN‘2 1 C DEL^VA )7 rG$*DTT) <"KTD*TFS1 

IF10FULT..CI) AoMgN 2 TO II 

inoci.i t *.nn iwu = i 

" 2 CM L AprTnUJMIUUlrAl) * “ ’ * * 

IF ( INC. EC.? .ANu. k T UGO .LE .R I GLM) IND3*U 
IF < TM~.C0.U1 GU 10 225 

IF 1 OTT .1 T .1 RO. .AND. IN0o.E0*1.) IND"4 ~~ 

1H IXC.ro. 41 AooiuN 4 JU 11 
I r ( t N C . r C .4 ) Gu TO 

P rt f L = AT AN2 ( (DbL^VA) t (GS*DT ) ) ^RTD^TBST” ’ * 

IF! ARSCFCLl UbURU) kULL=S I GN ( RL 1 1 ROLL ) 

IF(X1.LT.C.) I NO -3 

IFIPTrGP.LC.IUI^VAJ) 1 N0= 3 

I F ( rxr.EO.3) AGMGim J IU II 
IF ( TXC.EC.3) Go TO 3 

— 225 -CHl = MCGC/irD 

IFUMCO.FC.l.) VA1=VA*GGRT ( GS/ (CRL i#AMAG ) ) 

HHUH-0H1 

""V L » V A UrX^I-OHl/ 12 • ) ) 1 " 

VU"VUFCT“2 

pnr>rF=-vA/i.no 

— '225 V?=VAl*FXnMH-hF)/l2.*HS)) 

V22 = VAt rt FX°(-t H~MF )/{2**HSD) 

224 RP=V1 1**?/(GS*I ANRii*2U902QU0.) 

p*o -v 2?**?/ <Gj<-1mJKU«2 0902000.) 

1115 1=1.3 

1 115 T E n T )-UPT3HH'hR»%«Ukri(I)*lLr 

CALL IWlTVtTEf^ut.AlfXl) * " 

IF( I NC3 • P C • 1 • ) Go Tj 1114 
CALL ACCTElUPTl.UNI UCTH) 

mCTUGT.U) UH=1. 

STH--SORT { 1 .-C»mv*2 J 
CA1 L AnPTPUJNu,UFf AK*> 

jc< APC.LT.C.) 6 1 H=- 6 TH -- 

DO 1113 I = 1 1 3 

1113 TFM I) = L'PniT)M { RRk -KK* C TH )*URT 1 ( 1 )+Rk*S TH*UF< I> ) *ILT 
■ — CALL L\‘! TVMFrt.UKf*: .XU XI) 



10 :> 


1114 



FOPTPAN IV G LEVEL 21 


LNDSTR 


DATE = 73070 


li>/5:>/*4 


r 




1 

Jn* _ 
N> 

l 


TJX2S- 

0179 

0130“ 

0131 

0132 
0133“ 

0134 

0135 

0136 - 

0137 

0138 

0 1 39 — 

0140 

0141 
0142- 

0143 

0144 
‘014$“ 

0146 

0147 
0148- 
0149 
0X50 
0151- 

0152 

0153 
0154- 

0155 

0156 
0157- 
0153 
0159 
0 loO- 
0161 
0162 
-016 V 


0164 
C 1 65 

0166 

0167 

0168 

016*. 

0170 
' 0171 
0172- 


0173 

0174 
0175, 
0176 


CTHJtrrn^S UT(RIU*GS*HS*FTR* ( X ./VZ**2-r. /VI**:zr*RTD *FCTRT 

+ KCl*tCOS T NO<KL U* <HHt-HH*RTD/l2.*HS) 

I F C INC* ME *3) L)TH = L)Th~OTbIAS 

OAZ=t AZ^-AZIMJ^ILT 

DO 222 I H = l » 10 

IF< ( n-»360.*DA/J .GT.OTH) GO 70 2222 

“222 — CONTI MJfV ’ ‘ “ ‘ 

2222 DTH1 - IT H- 1 ) *3oO. *DAZ 

T F C (OTH-CThl-lt>Q.) . GT .0 . } DTHl^DTH l + 3o0 ♦ 

__ IF | INC.NE.3..ANU.IHX.NE.0. ) 0 7H1» IHX*0bO .+OAZ~ 

I F ( IKC.EO.l) GJ TO 2227 
CTO! F=DTFC“CTHl 

1 F { < 4 ESI CTniF)) .Gl .loO. > DTHi«OTHl+SIGN(3bO. , DTD IFF 

2227 DTHT= H THl 

I F ( IMUM:.?) UT Hi -= DT HI 

— £°0-=APS( RTt-l*SINOtKLU/OTH* * - - 

IF <APG.GT.U1 AKGU. 

PL2 = AQSI\'tAPG)*KJD 

PL2^PLl+KTUP^^tRL2-RLl) — 

IF< D L2.0T . ° L 3 ) RL2-RL3 
IF(°L2-LT . °L2rlN} RL2=*L2MN 

IF U NCI?. EC. 1. ) rANRU*$IND<Kl2)/CGSDl*L2) * 

l F( I NC 12 . EG.1 • > CKLL=L0SUIRL2) 

IF{ TNC12.EC.1. J FCrKX"CUSUtRLn/C0S0<RL2J 

on n 18 1=1,3 

1X10 TEPt I )=U* t ! )+kr>ULl 1 J 

2220 CA( L L N I TV(TFNrUt,l»Xl,Xl) 

CALL AXP < IJC ,nc l , I E«4 J - — ‘ “ 

ft3=0A\r v (ui» 9 uCJ*»F 

DPL=“KCN* ( &Vmv#;>F) 

CALL -LNITV<rrM f 0C2TAXtXl) 

CALL AO^TP<ULtUu^,Ll bSTJ 
OTL? = CAN’G (LC?UC1 J 

!F< INCn^F.l.) GO TO 1 X X 7 * - — 

I F l CT EST . EC *0. J GO TO 1X17 
0°L =— KfN* DEI 2*61 GNU. ,CT66T) 

IF<C7 r ST .IT .LM2 1 GO TO 1U7- • 

IFtOTFUtT.CThv.) oJ iO 1117 
CALL AHGrBUJt ,JF ,At) 

IFtABStXl ).LT.LiS3) bO (U 1U7- - • - - — 

IF(0E12.GT .CFLLlM) RLz: = 0. 

1117 IF< < CTF-OTFLW) .L f .IdU. ) GO TO 4441 

IT< nTH-LT.360-.ANL'. 1 0 TH-OT lU-OT HLM 1 ) . GT-.0. >— 60--T-0-6 

4441 . I F ( 1 MVFG.3) kC/LL= (RL2*UtM_)*lll 

I F C A3S(Pn.t ).ul .RL3) *0lL=SlbN<RL3fkniL> 

ALPH= ALF7 tKONPH RLH KEI — ROT ) . 

IF< TNC22.FC.1 .) ALPH = ALPH*CQSRLD/COi>Dl ROLL) 

COS^LC=CCSC(ROLU 

IF < INC.EQ-3) -ALPH=AIPH + KBLND* ( l.-DTHl/DTHl 1 *( ALFC-ALF7 ) 

CO TO 444 



f m t 1 ' £ x TV 'G LtVEL 


X5/55M4 


PAGE 0005 


to 

CO 

I 


1M”‘.IT>F) lNU-4 
IFtOTK. LT-DTniv) l,ML>=^ 

I F ( *NC. r 0.4) Abbl^M h TO II 
IF( INC.E0.4) bJ TO ***t»to 
HHI = w 

VA1=VA*S0&T (Gb/ { GKLI *AMAG) > 

V 1 - V A J ’ 

VI l = Vfi 1 

P nrr>rr^-VA / (LUU^CKH) 

'* G n rn 223 ■ 

CAtL A pf T H t U° t JK l X ♦ A kl> i 
on-A°r^<;F 

■ ' o^ = (vA^>-‘2ycrs>«-rA < 4RL>))*2.*rsiN(OEU/2: - n^? 

OLt. A or VM V * UR 1 i f V L A 1 J 

03 = Sir\(n*i t viAl ) 

* C2-3C4D3 

IFU A0Mn2M.Lr.U2L I M) INU^5. 

LRrc*mPGr-RlAbl*oOdO. >/H 

“ IFCINC4.cC. 5*) GO TO 41 " 

r n i ucelic* tu-u:)Cu-iuo>/ 2 .) 

innEi 3.1. T.C.l uuLJ=0. 

' IMnei3.GT.tPI/2«) J UtL3=PI/2. "" " “ 

C0 C L=-SI C,\[rr L } t oo> 
r FL = DEL*T f ST^DUtl 

TEST = l. 

POU®ATAN2C <ntL*v'A) , (Gb*OTM *RTO*TEST 
IF{ KOC.NF.O ..aNO. t.VO*. Eg. 5.) R0lL tt MJD*D0 *KVV *VL AT 
” IF( (ABSIRCLUKGT.RlbJ K0LL=SlGN(RL5tP0LU 
JF( T IF. FQ.K. AM). VA.U ♦ VLG) LG=i* 

LCCC- <nrG r -BIASjL*oUdO.)/H 

TF[LPCC.l.e.C.) luGOLOO 

POTPEF=-VA/LnO 
I f= C I«riO.FC.l«J GO 10 4Vol 

‘ WTT-HF2M*CT**2-R0?2**2)/U**GF2r ™ ‘ 

m M • l F.HTT) GO To ‘♦toi 
IF(H.GT.HSP) GO TO 

- = V-0 + '< SHM^h ‘ * ' 

= VA-vn^b) fbGL+KbblMLUO-lOCC) 

IM SB.GT. ICO. ) V3=100. 

* -IFtSB.I T.C.) bb= 0 * 

ALPH=ALFC + KALFM LuOCH.00 >/DLDDA*KOMPl*(RDTAEF-ROT ) 
ALPH = ALPH+KBlAbrt<*{ VI A-V AX) 

- !MT*c7.nC.U> ALPh- ALPH/CObD ( RULL ) ” 

IF( ^C3.FC.i..*NJ.*uMRULL) .Gf -RULLM ALPH=ALFF 
!FUN'Cn*FC.U.A^O. lHO‘t.fMt.5. ) ALPH=ALFP 

ST vitt /5TE ALT'VilUt IN ALPHA " " 

DAlF*ALPM-AlFULO 

‘ IFtABSIHAt cj.GI.OAlFLMJ ALPH- ALFULO+SIGN ( DAL FLM tDALF) 

•- Al F f nr = ALrF ' * * “ — - * 

IP t At FF.GT.At FNX) aLPH^ALFMX 


12 9- 


r TP II' AN IV G UVFt 


LNUilK 


DA IE <* /30/U 


PAGE 000 0 


“IFm PF.lf.ftlF™ ) ALPH=ALFMN 
CALL fRrTPlUNU.vIrtl 1 ,AkG) 

IFIAPr.LT.O.) k0LL=0. 

~ INI = I AD 
HF1=HF 

AZ M 1 = A Z I M 

Sfl'>=SP 

IFITNC.0T.2) Gu TO AA 03 

“CALL LAI O\GIUkT2,UA,UeO,RTDiLGN0,LAl T,LUNGT) 

GO T" A A 6 

CALL LALOKGIUC.UZf OEO |RTD»L JNO f LATT » LONGT ) 

— RFTl 1 ° K ■ — - — ’ 

IF ( M . L r • H c 2 • AND ♦ KU T MN • N fc -0 • I kDT2=-VA*H/ t,\TGGO-3I ASFU 
IF(RDl?.LT. c nTKNI ki»r2=kUTMAl 
-- JF(RCT2.GT .RDIriXJ rtUI2=ROTMX •• 
ic( H-LF.HTT2) FDT2=R0U 
1F(H.LF.HTT2) GF2=GF3 

~GFl’{-«CT+=nT2J/DT 

IFIGFl.GT.CF2) Gri=t,F2 
A1.PM=AL°H*( GF1+GSJ/AHAU 

- I W D 1 0 = 1 . • • - ■ ' 

GO TO AA62 

W°ITE(6,555) 

_FOR m ATI 1 - LANDING- SlIfi-BbYGNL) - VEH ICLE-CAPAd 1L I TY 1 1 

ASSICN A TO II 
1 * 0 = A 

-GO TO A A At 

IF I I OTHl-GSASlNl) (KL2 I*[)T#RTD/ (C0SDIRL2 ) *VA) ) »GT. 1GO . ) GO TO AAAI 
DELl=CTHl/RTD-Pl 

—-ROLL = AT AN? < IDEL1*VA )-r<GS*DT ) ) *RTl> 

IF(A0S(°i:LI I.GI.kLI) AUl.L=SIGNIRLl»ROLL> 

IF I I ADA . EC. 3. I GO 10 6o 

- - HI =H 

I FI INC. 60.1 ) INUA=3. 

GO TP AAA 

- HF=hF+H 

h 1 r H 

GO TC AAA 

DP AAA7 1=1,3 - 

TGM I n*UPTU) 

CALL LimVITEM, 11*12, XlrXl) 

ALFPLC=ALPH 
GO TO A A A P 

ftLPH = AI. i'CI+KHM FI >-''(« UfRFI-kOT) - — - — - 

IF! VI X.GT .VI M2) GO TU 
ALPH=ALFCl-KTR*t VLM2 -V1X) 

AL PH= ALP H vKCMFT® l kl) ) kF T-kOT ) ; 

JF< ALPH.LT.ALF7J AL FH=ALF7 



130 


. >: jr/ss tv C trvFl ?t LNOSTR 

si*,? lcfVIX.l F'.VtMjJ GJ rO >‘ j 5 o 

l = (tsrr. ,: C.'3*J T£ hi 1 = 1 bi r 

• IFt'IFA.GE.CFL IK) Tt;.>U=rfcST 
'* p yn.-q\= { v 1 yx+z-Vf** l ) / (2* M ACTR) 

A 3 C» -R 7 rr *P / n TPGGn 

TfiA n c.rr,K) akG=i. 

pl 7 = i p p r,> *r id 

ir<cj 7.1T-W17MIV) KL/-RL/KN 

I F ( ° 1 7 • GT • D IT1) M.7-KLT 1 

Rr«tL = Rl 7* T T ST I 
G n TO 4 A 4 

-*>555 0 "“~S*$ T V r (°l d>/ (C0SDIRL5) *VITO 

I F ( nM « ; iT ) UM=UMLM 

IPUrMfrpT ) .CT.utU u«“UEL/DT 
• — - - F^t L “ A N2 ( ITMvV 1 X ) »G:> ) *TEST^RTD 
GO JP AAA 
FMO 


l 


'T7*> 

0279 

0230 

0281 

*“• 0232"" 

028? 

02 AA 
0 2R9 
02 Ro 
028 7 

02R8 

C? 89 

0290 

0291 

0292 

0293 


DATE * 7 3078 


15/5S>/4A 


PAGE 0007 
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IMT 


DATE = 72X58 


10/55/38 


PAGE OOOX 


U001 


SUbROUT INE IN IT 

t J02 


IMPLICIT *EAL U-N) 

Q0o3 


lNTEULR l ,M , J 1 1 ASStC 

OJO^ 


01 HENS ION k(3)*V(3),UZ(3)»UR( 3) » UNI ( 3) »UV( 3) t A t 165) » B < 51 

0005 


DIMENSION Ufc(3)»URZ(3),LIRT(3) 

U 0 06 


EOUl VALENCE (R,A) ,( V,A(4) > ,(UL,A(7) ) , (UZT, A (10) ) 

ooor 


EOU IV ALENCE l UR, A( 1 1 ) ) ♦ ( URT , A( 14) ) , l URTT , A ( 17) > , ( I LON , A ( 1 8) ) 

occa 


DATA YES/4HYES / 

oOOS 


tf(M.Ej.l) CO TO 11 

0J10 


M=1 

OOll 


Rl=2052o043. 

Oj12 


RP=2C85587 3. 

0013 


ATK=1 »/60 • 

0014 


KH=1.-(R1/KP)**2 

0015 

11 

WR IT£< 6, 21 ) 

Ou 1 6 

21 

FORMAT ( • NEm ENTRY FILE? '1 

0017 


RE AD ( b » 22) C 

00 Id 

22 

FURMATl A3) 

CCL9 


II- (C. EC. YES) GC TO 12 

00 2 G 


DEFINE FILE 7(33, 8C,E, IASS) 

0021 


RE AC (7* 1,60) A 

0022 

12 

hRl TE (6,600) 

0023 

600 

FUNMAT ( • INPUT:') 

0024 


READ (5 , 100 1 I 

0025 

100 

FORMAT (14) 

0026 


tr ( I . cO . 0 > GO TC 2 

0027 


IKI-GE.O) GO TO 101 

0028 


RE AO ( 5 ,200) 6(1) 

0029 

200 

F0KMAT(G10.2) 

0030 


GO TO 5 

0031 

101 

RE AO ( 5 , 20 0) (B(J),J*1,5) 

0032 


00 3 J= 1, 5 

0033 

3 

At je i -i ) =b ( j) 

0034 


GO TU 12 

0035 

5 

A(-I)=bl 1) 

0036 


GO TO 12 

0037 

2 

CALL UNI TV(R»LR«RL»RLSQ) 

0038 


IF (KL.GT.RP) GO TO 7 

0039 


IF (JZT.EC.O.) R(2)=U2(11 

U04U 


R11=K( i)+Rl/SCRT( l.-KH*SIND(R<2) >**2 1 

0041 


R ( 1 ) = Rl l 

0042 


R( 2) =C. 

0043 


R( 3 )=0 • 

0044 

7 

IF(V<1) .LT. 90000.) GO TO 8 

0045 


V< 1)=V( 2)*S1N0( V(3) ) 

0046 


V(2)=V(2)*C0S0(V(3) ) 

0047 


V(3) =0. 

0048 

8 

IF (UZT.EQ. 1.) GC TO 4 

0049 


CALL AX 8 ( V, R, 7 EMP ) 

0050 


CALL 0NITV(TEFP,LNI,XI,X1) 
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IN IT 


DATE = 72158 


10/55/38 


PACE 0002 


0051 


CALL AXBIUR.UM.UV) 

0J52 


CL=COSU<UZ( n » ■ 

0053 


$L = S t ND (UZ ( 1 ) > 

00 54 


CAZ=C(jSU(UZ<2) > 

00 55 


SAZ= S 1ND( UZ l 2 ) ) 

U05o 


uzm=cL*(uvu >*caz-uni< d+sazi+uri u*sl 

00 5 7 


0Z(2)=CL*(UV(2 >*CAZ-UNI<Z}*SAZ)+UR(2>*SL 

00*0 


UZ< 3 )=Cl*(UV( 3)*CAZ-UNI{3)*SAZ)+UA(3)*SL 

OOdS 


UZT = 1 . 

OOtC 

4 

IHUTI.Eg.1. ) GC TC 31 

0061 


IFIUKTT • E 0 . 2 • ) GO TO 40 

0062 


CALL AXtKV.R.TEMP) 

0063 


CALL UNI TVI TE PP ,UNI , XI , XI ) 

OUo’t 


CALL AXU lUk. UNI * UV ) 

0065 


LUk=CCSU!URT ( l )/ ATK) 

O066 


SDK = SI NO l UK 7 I 1) /ATK) 

OOo 7 


CCR=CUS D (URT { 2 )/ ATK ) 

0068 


SCK=SINU(URT<2)/ATK> 

0069 

, 

DO 32 1=1.3 

0070 

32 

URT ( 1 ) = UR< I )*CCR ■*SCR*(UN 1(1 ) *SCR+UVUI *CCR ) . 

00 7 1 


0RTT=2. 

0072 


GO TU 40 

0073 

31 

CALL AXB (UZ. UR.T EMP ) 

0074 


CALL UNITV< TEMP.LE »Xl ,X1) 

0075 


CALL AXB( UE.UZ.URZ) 

0076 


CLAT=COSO(URT (21) 

00 7 7 


SLAT=SINU( URT<2>> 

0078 


CLUN=COSO( URT l 3 )•* I LUN 1 

00 is 


SL0N=SIND(UKT(3)-1LCN> 

0080 


00 33 1=1,3 

OUtil 

33 

URTI1 )=UZ( I)<'SLAT+CLAT*(URZm*CLON+UEm*SLON> 

0082 


URT T =2 . 

0083 

40 

WRITE (7» 1,6) A 

0084 

60 

FORMAT (5E14.7) 

0065 

6 

FORMAT ( 1 P5E14.7 ) 

0086 


RETURN 

0087 


END 
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MAIN 


DATE = 7306/ 


11/41/37 


00C2 

0003 

.. . * 


♦ 

0004 

nnns 


0006 

* 


. „* 


♦ 

0007 

* 


* 

0008 

nno9 


00 LO 

* 


* 


oou 

00L2 

_.0013_ 

0014 

0015 
_ 00 16 

0017 

0018 


INTFGCW T < 1 1 . DC , I n 

D I !•' E M S T ON Will ,« ( 3J , V13 J ,VA(3) ,Jt\<3) ,Dl 3) .LI 3) ,G4 (3 ) , 

TFI'P l 31 *Tfcf'Pll3I *UL 1 31 ,LIVA( 3J ,U3 ( 3 1 « UZ <3 ) *M 3 ) »Z l 81 t 
CZ ( 8 ) « A ( 1 6 5 ) < UkOI 3 1 , UCk ( 3 ) , Ut 1 ( 3 ) . Ut 0( 3) 

DIMENSION CCttSJ ,X(o) ,IHRUSTI.>J 

EQUIVALENCE I UZ. A( 7 1 1 , I UtiTP , A 119 ) 1 

EQUIVALENCE IPClisuAlcU). I UMAX* A 1 22 ))» I HMIN , A 1 23 )) * 

t ,r, LL »A (24) ) t ( nOCUA, A l 26) 1 » t LL J* A ( 2o )) * ( i)T 1» A t 2 / ) ) t (LONO*A( 18) ) , 

( T. M AX t A( 29)) ,ti«iuAJUr At jJ} ) , Ut-iTk , A< 40) ) » IA131) tNO^LT ) , . . 

(A(32)»DT2)«(A(j3J» HMN<: ) , ( a t 3 A > , PC I K2) 

EQUIVALENCE <ZfVJ*lZl4)th) , (CL A, A ( 15*») ) t (COA , A( 15 5 ) ) » 

(AL p H»A(156))#(l KAo iAtl67))t(i v U)Ur.All68J) i IM . A( 159) J * . . ... 

(KW IF *A ( 163 ) ) 

DIMENSION HM7I ,V«1 17) ,UN1I3) , VW13) 

DI. V EMS I CM E112J 

EQUIVALENCE <AIa7),.iF),IAIao), AZu ) , t A I 3tt ) »KC LA ) , I A l 39) >KCDA) , 
IA(56),TW),(A(57),TA), ( AU>3) f VGS T ) , ( AI54 ) ♦mROD ♦ 

( AI 55 ) t KApr ) . - 

DATA HW/30.,60.,23.»20.»14..10.»0./ 

DATA Vrfl/150., l30.»40.*40.,97.,97. t 28./ 

.. OATA CQ/S+O./ . - . 

DATA X/8*0./ 

DATA Z. 02/ 16*0./ 

DATA LG»SB/0.,0./... 

I H = 7 

IFIM.CQ.l ) GO TO 1 


U) 

to 

T~ 


.0019.. 
0020 
0021 
>022. 
09 1 3 


GS=32 .2 . .. 

MU= 1 . 407654F1O 
W I E=« 0000729211508 
KQ= 3. 049F-6 
J0= 1 o2 3 .46E-6 


CO 24 


P 1=20926043* 

0026 


RP = 208553 73 * . „„ . 

C026 


P F= 2092 5 738. 

0027 


KH=l .-{91/RP>**2 

0028 


APFF=7320- 

0029 


DC^l 

0030 


RT0=180*/2.141592o6 

00 31 _ 


ATK = 6C.**RTD 

0032 


Xl = 0. 

0033 


PCT--3. 

0034 


fnop=o. 

0036 


r=o. 

00 36 



00 37 


CALL IN IT 

C036 


PE AO (7’1,21) A 

0039 


DT = DT 1 

0040 


PCTR = PCTP l 

0041 


CAZ^=COS0(A2W) 
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MAIN 


DATE « 73067 


11/41/37 


004? 
Cv *3 
C ) .4 
0045 
00 V6 
00 V 7 
0048 
OOtv 
OO-'O 
C 03 I 
0062 
0053 
00 6 A 
00 6 5 
0 0 '3 b 

0057 

0058 

0059 


"" OC 60 " 
CG6 t 

0062 ■*” 

0063 

0064 

' CO Ci 5 ■ 

0056 
0C67 
0068 "* 
0069 
OC 70 
00 71 
00 72 
00 7 3 

0074 

0075 

0076 

0077’ *“ 

007 8 
CO 79 
0080 
0081 
0082 

0*103 ” * 

0084 


SAZW=SINOUZW) _ , 

" $^T=siNOU pn n 

C^T-CO$D(V>pnr) 

D n 27 T = U3 
R ( I ) - A { I ) 

27 V ( I ) = A ( I + 3 ) 

21 FORMAT ( 5F 14 «7 J 

ce n u.»crsn<puu 

S^nileSl N’C(? n LD 
WIE = KMF*ME 

CALL LNI TViRtUROtAl fXl) 

CALI VFCSCl (UZiwict W ) 

♦ CALL AXA ( UZ tOKOt TtM P ) 

CALL L\HV(T c .<i»,UfeO»XltXiy' 

CAI L AXEM V.I^O.le.-IPJ 
CALL LMTVtTPHP, JCK.XUXl) 


51 

FORMA T ( l HC t 1 1 l Mb 

RANGE 

ALTITUDE 

VELOCITY 


*0T ACCEL.' 1 



... 

►— . - 

’ 1TE(6»52I 




52 

FORMAT (• 0 

LAT 

LONG 

LATRANGE 


* WACH') 


1 IFIDC.F0.2) M=rt«-.3*UT*M00r 

* IFinC.EQ.A) «=M+.5>»UT*MU0T 

CALL LNtTV(R,UK,uL»RL$y) 

• CALL AX0(VJ,»,1bMP) 

CALL VSUiJ ( V ?Tb V » 1 cMPl ) 

CALL AXML'Zt'JK.UMF) 

' ‘ ’• CALL LNlTV(TEMPrUtlfXlfXl) 

CALL AXt3(LT,Utl,0NH 
CALL AOrTIMU=,UZ,CI> 

H=FL-Fl/SC t, T < i.-Mi''Cr-i‘«‘2) I 

1F( INCF.EC.l.l H-ZIti) 

DO 1 0 l I I = l » I H 
' |=lH-ll+l 

IF((.30Afl*F).LT.(1000.*HW(I))) GO TO 102 


101 CONTI MIF . 

■■“'102‘ CL A") = ( . 0003U‘tfc*M-Hrf ( I J )/ ( HW ( I ) -Hrl( 1+11 ) 

VW2 = tVVU< t I+GkAOv (Vtil < i j-vwlt 1*11 n’MWf-/.30< t 8) 
if(t.ct.<tv»+ta); cj ro 203 _ 

IF<T.GE.TVO Vr(2=VGbr"H.-C0i0l(T-TW}*3o0./TA)l 

203 31 T=l.3 

31 V«( t) -VW2M (ti;»l( 1 )*CAZW+UEl< 1 }*SAZW)*CWRT+<J:t<n*;>WRT) 

CAM. VAO I V'a *Tc'3F 1 » V A J 

CAI L LNTTV(VA,UMA,VAL,VALiy) 


0005 

0086 

0087 

0^35 


UEs^E/^L 

* Bl=<l.M<l.-5.*LT**2>*J04'(.5 + (10.5*CT#«‘2-7.)*CT*#2)*K0* 

* UF’M'2 )»UF.*A2)w(-M01 / iRLiw) 

B2 = (2.*J0*(2.-i'*.«‘Cr**2/:>. ) *UE**2*K0 )nJE**2* (-MU) *CT/RISQ 
00 2 T=l,3 


RD 

ROLL 



FPPTPAN IV G LEVEL 21 


MAIN 


l)A IE '» 73067 


11/41/37 



c - 

END GRAVITY CAllOLAUUN 

0090 


CALL AT M 0 $ (H, KHU , Viu'JNO ) 

CCS 1 


MACH=.3043«VAL/VSUUN0 

0092 


Q=r<,5* ( P HQ *KBHJ)*VALbU 

0093 


CALL AXMUVA,UK, IbNP ) 

0094 


CALL LNITVlTEMP,uL,Xl>Xl) . . 

0095 


IF(OC.GT.l) GU 10 ill 

0096 


IFINCSTR .EO.O.) GO TO 111 

009 7 _ 


.CALL L'NITVIV. rLhP.VMAG.Xl) 

0G98 


CALL LNDSTRtP*V»UtO»VAL»VMAG,AZIM»H»f\JLL»ALPH,/ 
* B(3) , B(4) ,6{5),t>(o) » B { 7) , B { & ) ,0(9) » b ( 1 0) ,8(111 

0099 ..... 


_ 1F{°CT#F0-C • ANl>* GUTP*Nb*Q) . WKl Tb ( 0 , 26 1 b . 

oioo 


CALL AFP DU , A Kc F ,R , V , VA ,CL A ,L DA, ALPH ,kOLL ,XACH 

0101 


CLA = K Cl. A<«C L A 

. 0102 


COA=KCCA*CCA . .... 

0103 


SPOLL = S l NO ( POLL) 

0104 


CROLL=CQSD(P n LL) 

0105 


SAl PH-S I NO { ALPH 

0106 


CALPP^C^SDI ALPn) 

01C7 

111 

01=0*CDA/M 

0108 _ 


CALL VECSCt. IUVA ,-Ul • 0) . -. .. 

0109 


CALL AXG ( UL ,UVA, N3) 

0110 


00 3 1=1,3 

0111 

3 

Nm=K3m*CR0LL+0L(l)*SR0LL 

0112 


DLr‘0 = C*CLA/M 

1 0113 


CALI. VSCSCL(M,OLUO,L) 

I" 1 ' 0114 


DH 4 T= 1 , 3 

S 0115 
1 0116 

4 

TH^US T ( I ) = ( TMAG/K) # ( N{ 1 ) +SALPHHIVA C 1 )*CALPri) 

C 

IFIPCT.LT.O.) GO rj 121 

0117 


DZ ( l ) =LI l ) +D( l)+*THi<UST it 1 ) + GRl 1 ) 

0118 


0Z(2)=U2)*D!2H-ThK0i7t2) + GKl2) 

0119 


02(3) =L 13 )+0t J)rfnKUST (3)»GM3 J „ . 

0120 


"" ” DZ(4) =Z(1) 

0121 


0Z<5> =Z(2) 

0122 


DZ<6) = 2 1 3 > . * 

0123 


"iFlINCF.MF.l.J GJ TO 333 

0124 


CALI AOOT 2 t OZ i UR , Ah ) 

0125 _ 


0Z(7)=AH . , 

0126 


0Z<8)=Z<7) 1 

0127 

G 

333 

CALL CtFEQt8,DT,UC,DZ,Z,J,t)QrX)_ - .. 

0128 


" IF < DC .LT.9) CD Id 1 

0129 

303 

OC=l 

0130 


IFOCT.LT.C) Gu TO 121 . . .... 

0131 


PCT=PCT+1. 

0132 


IF( T.GE.TMAX) ENOis=l- 

0133 


IF(H.GE.HMAX) tNdK=l. 

0134 


IF ( H. L6.HM IN) cNd*=l. 


B< 2) t 
LG, SB) 
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MAIN 


DATE = 73067 


ll/Al/3? 


PAGE 000A 


0135 

0136 

- 


013 7 



0133 



0139 

0190 

0191 



0152 
0 1 h 3 

121 


0194 

0195 

*■-* 


0196 

0197 


. 

' 146 * 



0199 



0130 

0151 


-* 

0152 

0133 



ri:>9 



0155 



0156 

”0157“ 



0153 
0159 
’ C 1 60 
0161 

■ - 

-- 

^ C162 

_ „ 

* 

0163 

501 


0 lo4 

‘"C165’“ 

13 

- ~ 

0166 

0167 

0168 

12 


0169 


* 

“ 0170 

“25 " 


0171 

0172 

26 


“ 0173 

0179 
Cl 75 

11 



IFUNOF.FO.l. ) GO TO 121 . .. 

IFtH.LE.HWN2) INuF= l . 

IF(H.LE.Hw\2) Ot=OtA 
IF(H.l F.H'.*N2) PCTK=FLTR2 
IF t H. LE.H. M N2 ) PLT=PGTR 
I F ( H . t s .•t‘ / N2) A t d ) =H 

I F ( H.LE »H’lf»2) CALL ADOTB ( V, UK, Zl 7 ) ) 

CAI L VAC( t *0»TtfriP) 

CALL VAGI r r,, P»lHA06T , ItMPl) 

CALL LNTTVI TKHP1 ,TEMP,AMAG,X1) . 

CALL AnfTRtU- ,U n.O,AKG> 

RANGF=AT**ASCG!>IaR.>) 

CALL AC°T R ( 'J e «UC A » aHG) . 

L A T 3 f i C = A T K *. A" 1 u 

CALL UNI TV I V* TEMP » Vfc L » XI ) 

C A l L ADOTRI V, Jh,tLOUT ) 

LAT=A B STN<CT)*r. 10 
CALL AXBIUZfllAfl EMP) 

CALL LNITVITFMPtUbl tXliXl) . 

CALI. AXP(UPl.JLU,IcKP) 

CALL LM T V ( TF,ii*» TLMP 1» AftGt XI) 

IFtXl.EQ.O.) Ar%y - 0. . 

[)t.ON= APS I M A B G) 

CALL ADOTP (UZ * ) £MP t ARG) 

IFIARG.GT.C.) OLON=-DlUN _ 

LCNO = LCNO + lR' UiN-wleMlPfiTD 
IFCinPt T.EC.l.J gO TO Id 

VvPT TE < 8» 501 1 T f n. ANoe r ri t VAL»KOOT , AM AG , Q , LAT t LONG* l AT A NG* ROLL* 
MACH,e 

FOPMATI 1P6E12.A) 

IFtPNOR.EQ.l.) GO TO 12 _ _ _ 

IFI °C T.EC .PCTR ) GO TO 12 
IFtPCT.LT .C) GO Tu 12 

GO m 1 

PCT=0. 

W>JT£<6*25) T, RANGE *H *VAL*RDOT *AMAG*U «LAT *LONG*L ATRNG*KOLL * 
MACH 

FOPwATI/lPfcFU.H/ot 12. A) 

Fm MAT( IP6E12.<*/6E12.<*) 

IF (ENOR.EC.l.) GO TO 11 

G" TO l 

STOP 

END 



-p u 




9 


7 


rtEAOY 


...n r v< * 1/ 'U-uu* ^ * 1,1 

J Ul=* M I f* ) +»■*!'? A M 11 1*. 1 l f +1»K )-LU(l » ,J s ) ) 

- iV -/ :(I ,. 4 +j) » . - ^ o * c ^ ‘ < 1 +3t:; t< +D* 

Ci,»,: •- :* 1 ) 

) T tl + ■>*•% '.S* (\. - J ) 

ir < ' 1 1 r ‘ J * »■!'•) 

CUl-O'.TE : ^ 

4 > <r /••* 

r u * t i«i.r • . l ) i 5.> t * » 

Ir ( ( , M ,♦ o T » / ) * *0* i- 7 • 

II Mi 

r< £n> - xr( i) »» {•vjh-d+u'UU) 

cn i-(- i >- : /Ok. ^MM1+IX>*C JO+AKkF 
j r *y J«M. 

L»'H/u.-J 
00 •' •• = ! t *> 

■< -•>- +1 

in*ir. v-.^o^rtK )> oo to 9 

conn me . , 

3 ,i a )1=(-0H-*«''IMK >> >/(l«0HTU +X >-U0HT < K .) > 

i)CLl ♦ n :;i 1U *J ) A MU>l*tOCLl(K + 1 • J :>-*JCLl <K '0 > 1 

*)CL l V— :V,LU< tj '+1 )«j*AUI*(JCUUK*+1>0 +1>-0CL1<* 

UCL= *'-t *'t -CL1 ii-.'Oul ) 

cua= ;> mHJ % ** v>**> 

CLu-Ct A+U«‘.» 

CUAsCOirMSJC 3*S i+L'S^C!) 1 0) «AK£F 

•<t£f I \M 

£N') 


> J :+l) ) 


u> 

03 

1 



> '.I 


20 


ri 

2J4 


2 . 
2d 


20" 


V 


30 

31 


w 

VD 

I 


c 


1 

2 


a 


o 


3 


S'M’OTril •• Tf.) X) 

: > ■ > * 0/( 0 • *t >#>.,*< i> ^ / C 4> 

■ )*> 4 J=l P‘v 
P ) T C eU # "'I tj f. 3) > mA 
CGr.b 
Cl = 1 

c j=#b*n r 

v* ) X ' 20 . 

COS. 2V 'J liiJ' 3 
C 1=00 

c-»=. o*ur 
v- > r« • <o 1 

C<!=1 • /i) / into /ttJ * 7 
0) r . 214 
C*>=# J *> 

C i - • 3 
C2=.b+OT 

<c j)=co*n'<j) «0T-ci* , )( j) 

0(JJ = XJ) +3'*X(J) - C2*UZ(J) 

7 ( J) - ( J> + X ( J) 

G ) 1 1 • (30*»i ) ** l A 

f-=- +H. /2 

V=* 41 

»■< ; TO K i 
i£"JU 

?* Ji i*?00T I Nh A ) 'AO** (AH.IOrV^J' J»NU ) 

: 3 l P( 14 ) J 4 )ri 1 .|«)U 5 J pM*>H 0 ( 14 ) 

£).\ r ' U 00/1 *btib r 1 *nc i • 2 Kb * 1 • r.bf 1 «G:.br V- t- 4 f 7 .Oii 4 * 

* b * lt£4 • b ♦ 2 *' t . /- ' * 3 * 2 .4 1 2 ♦ •„** ♦ 1 • i 1 . * ♦ 0 • / 

U ‘ f *- «H/1 (J.mrV ''vOt,} U« vJ>2‘ !)■' b>2t0*bt*>r 1G0 *bb # 

* ltL)»'>5r2> •< •*>* 27'J H#bbr21 b*210* b>2‘3 .lb/ 

u \ r ■ ( t')/ . 0 lb f « 11 2 p * 0 JL t * b b f » b 3 » ' • “ “ • b 4 r — ♦ U ^ r J • f • u . 2tJ 

♦ *0 *1 m)« f "* (i bb/ >■ 

OU' 0 tr* 0/1 • • 1 • tt3 , ^E—9 r 2 . 43(n.-P f y •ti2V‘£'"‘J • 4*97*’ £-/ • 

♦ 3*17«rh-Of2.U 1 t-b » 2 ..:>10yc>4 r 7 . b943*-* W 1 4427bh-3 r 1 .32 ibt-2 * * 
t »:m Gob*. *6 v2#1.2'b/ P 

S3 iT '*_/'v3.> *0 ./ 
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APPENDIX E 


NAV-I GAT ION F I-LT-ER CHARACTERISTICS 


The navigation filter of Ref. 112] was used to generate 
the precision DME results. This is shown in Table E-l . 

The Triangular Formulation of a Square Root Kalman Filter 
[16] was used to generate the low precision DME results. 

This is illustrated in Table E-2^ A dramatic increase in 
computer speed (nearly 4 to 1) resulted in using the Square 
Root formulation. A nearly 2 to 1 increase in precision 
characteristic of the square root formulation also resulted. 

A thirteen state filter was used. Two states were 
added to the filter of ref. [12] to estimate the 90 m bias 
of the low precision DME's. These states are illustrated 
in Table E-3 . The corresponding state transition matrix 
is in Table E-4 , and the driving noise elements are in Table 
E-5. 


The landing aids are located in Fig. E.l. The two 
DMEs are rather close to the runway and probably should 
be moved out. But this baseline was chosen to provide 
results comparable to those of ref. [14] . 

The DME model, Table E-6, is from both ref, 112] and [14], 
The 90 m bias required the addition of two filter states to 
estimate this bias. Without the estimate, divergence 
resulted and even with the positive definite guarantee 
of the Joseph formulation, which was used in early tests, 
the covariance matrix became negative] 

The barometric altimeter model, Table E-7 , is from 
ref. [14] . No bias proportional to altitude is present 
(scale factor error) and probably should be corrected. 

The radar altimeter model is in Table E-8, and 
the ILS model is in Table E-9. 

As in all these studies, correct modelling of the 
navigation is vital. This is true both in the filter 
equations and in the error driving model. For example, 
the filter will drive the altitude close to the .61 m 
value of the radar altimeter near landing. If terrain 
effects are greater than this value, the simulation will be 
over optimistic in the altitude performance near landing. 

Finally, the IMU characteristics from ref. [13] as 
summarized in Table E-10. The 2/3 meru gyro in this unit is 
a high quality unit placing this system near the front in 
the state of the art. 
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Table E-l: Kalman Estimator (Joseph Form) 

Used in Precision DME Work 
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S T s = P (S = upper triangular root of P) 
f = S T h 


a = r + f f 


S + = S[I-f f T /c0 1/2 


where [I-f f T /a] 1//2 is upper triangular root 

of I ] , derived by Cholesky 
, decomposition. 


x + = x + Sf Az/a 


Table E-2: Triangular Formulation of Square Root 

Kalman Estimator Used For Low Precision 
’ Navigation Study 
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— — 

Variable 

— - , ... ■ * — 1 

Sign Convention 

X 1 

Error in east position 

Positive if indicated position 
is east of actual. 

X 9 

Error in north position 

Positive if indicated position 


is north of actual. 

x 3 

Error in altitude 

Positive if INS indicated alti- 
tude is above actual. 

X 4 

Error in east velocity 

Positive if indicated east 
velocity exceeds actual. 

X 5 

Error in north velocity 

Positive if indicated north 
velocity exceeds actual. 

X. 

Error in altitude rate 

Positive if indicated up 

6 


velocity exceeds actual. 

x 7 

Platform tip about east 

Positive if platform is rotated 

/ 

axis 

positive about the east axis. 

X 8 

Platform tip about north 

Positive if platform is rotated 

axis 

positive about the north axis. 

X 9 

Platform azimuth error 

Positive if platform is rotated 
positive about the up axis . 

x i ^ 

. i 

Vertical acceleration 

Positive if it induces a 

10 

error 

positive altitude-rate error. 

X 11 

Altimeter error 

Positive if measured altitude 
exceeds actual. 

- X 12 

DMEl Bias 

Positive if bias exceeds actual 

* V 

hi 

DME2 Bias 

Positive if bias exceeds actual 


* New 


Table E-3: State Variables Estimated by the Kalman Filter 
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rr T 


*7,7 - 4 ' 


*8,8 lr 


$ » I # 

*9,9 ' 


*4,8 " ' 4v z 
*4.9 = 4v n 


*5,7 " 4v z 

* 5t9 “ - Av e 

* 6 ’ 3 = 2(g/R)T 

*6^7 = - 4v n 

*6,8 = 4v s 

* 6,10 _ T 

4 > ? 5 - - ( 1 / 11)1 

* 7 ; 9 - - 4e n 
*7,3 = 4 V R 
*8,4 - (1 / R)T 
*8,9 - 49 e 
*8,3 " - 4 V R 
*9 4 “ (tan L/R)T 

*9^8 = - 40 a 
*9,7 40 n 


o = ” C v tan L/R )T 
t o e 


10,10 = 1 - 
’ll, 11 = 1 - 


(v/d gz )T 
(v/d l5 j.t) T 


~ 1 - T/400 

** ' = 1 - T/400 

9 13.13 V 

* New 


TABLE E-4 : NON-ZERO ELEMENTS OF THE STATE TRANSITION 
MATRIX 
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°4,4 ~ 

|Av e | v 

°ASF + 1 

;4v n | 

v a lo 

°5,5 " 

|Av n i v 

° ASF + 1 

l iv e l 

* 

°6,6 

2 > Av horl 

l V t 20 A6 + 

(o ABIAS /3) 


^ ^ 2 

^7,7 “ At T B (°GBIAS +a ADSRA g ^ + V °ADIA 

+ l A0 e l A6 total a GSF + At(2v/d 6n )a 6n 
Q 8,8 ” At T B^ a GBXAS + °ADSRA 9 ^ + ^ Av n^ V °ADIA 
+ l A6 n! A0 total °GSF + At ( 2v / d 6o )a 6e 
°9,9 = At T B (a GBIAS + °ADIA g2) + ' AV J V °ADSRA 
Q 10,10 = At(2V//d gz )a gz 

Q llfll = At(2v/d alt )a 2 lt + |Ah|h s c 2 emp + |A(v 2 ) 

* Q .. = 2 GBIASD 2 AT/400 = Q 

* °13,13 ~ Q I2 t 1 2 

* New 


2 2 
v o 
s sp 


Table E-5: Non-Zero Elements of the Noise Covariance 

Matrix 
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APPROACH 

DIRECTION 


DME 1 
©-<— 


Runway 


Localizer 
Station (4 ) 


3 Km. 


V2 


■ 153m 




GLIDE SLOPE 
STATION (6) 


Fig. E.ls Landing Aid Locations 
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The measurement variance is assumed to be 112] 



90 bias is added to DMEl, measured range -90 is 
added to DME2 and a random error of 7m is also added. 

The filter estimates these biases in states 12 and 13 
assuming an exponentially correlated random variable 



Compensation for non-linear elongation of the measured 
range are made to the actual measurement and assumed 
variance as in ref. 112]. 


Table E-6; DME Model 
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The altimeter has both a bias and a correlated 
error. The bias variance is: 



No bias term proportion to altitude (scale factor) is 
presently added . 

The uncorrelated err.or is also added as: 


2 * 2 

0 2 

0 = K 3 h 

+ 3 

„ „-8 

K 3 - 

6.25 x 10 


This bias and uncorrelated error are added to the measurement. 

The filter assumes a measurement noise of the 
uncorrelated value. 


Table E-7: Barometric Altimeter Model Ref. [14] 
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The radar altimeter error is altitude dependent 
as shown below. 



At present this is added to the measurement as 
an uncorrelated error. 

The filter assumes a measurement noise of this 
value plus an additional 1 m, additive value. 


Table E-8 : Radar Altimeter Model Ref. [14] 
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The localizer (azimuth) error is range dependent 
as shovm below: 


range p (meters ) 

a azimuth (mr) 

0 < p < 1067 
1067 < p < 8330 
8330 5 P 

1.4 

1 , 4+4 .2 (p— 1067) /7 2 63 
5.6 


The range dependence of the glide slope is shown 
below: 


range P (m) 

a elevation (mr) 

0 5 p < 1067 
1067 ± P < 8330 
8330 5 P 

1.4 

1. 4+1.1 (p-1067)/7263 

2.5 


Uncorrelated random numbers of this magnitude are added to 
the measurement. 

The filter also assumes a measurement noise of this value. 
If measurement is not within 2 degrees of the runway azimuth 
or elevation angle it is rejected. 


Table E-9: ILS Model Ref. (143* 
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Source 

Units 

Value 

gyros 



Bias 

deg/hr 

.010 

g sensitive 

deg/hr/g 

.015 

IA Misalignment 

sec 

60 

Scale factor 

ppm 

50 

Accelerometer 



Bias 

yg 

60 

Scale Factor 

j 

ppm 

34 

XA Misalignment 

/ — • 
sec 

40 


Table E-10: Inertial Measurement Unit Model, 

Kearfott KT70 Ref. [11] 
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APPENDIX F 

Aerodynamic Characteristics for the 040a Vehicle 


MACH ALPHA - CL CD L/D S.BU) 


o.aohMo 

0.0 

0. 02750 

0.30000 

2.00000 

0.00162 

o.3unuy 

4.00000 

0.14075 

0.30090 

6.00000 

0.20550 

0.30009 

a.onoou 

0 . 26626 

0.30000 

10.00000 

0.32675 

0.30000 

12.00000 

0.39125 

0.30000 

10.00000 

0. 45562 

U. 30000 

16.00000 

0.52000 

0.30000 

1U. 00000 

0.50200 

0.30000 

20.00000 

0.6^400 

0.50000 

0.0 

0.01550 

O.bOUOO 

2.00000 

0.004-57 

0.50000 

4.00000 

0.15425 

' 0.50000 

b.uon no 

U.222O0 

O.bOUOO 

fl. ooooo 

0.29075 

O.bOUOO 

1 n.oooou 

0.35725 

0.50000 

12.00000 . 

0.42375 

0.50000 

14.00000 

0.49187 

0.50000 

16.00000 

0.56000 

0.50000 

i«. ooooo 

0 .62600 

O.bO'JvjO 

20 . 0000U 

0 .69200 

• 0.70000 

0.0 

0.01635 

0.70000 

2.00000 

0 .09U67 

0.70000 

4.00000 

0.16500 

0.70000 

6.00000 

0.25767 

0 .70000 

6.00000 

0.31033 

0.70000 

10.00000 

0.38150 

0.7000 0 

12.00000 

0.45267 

0.70000 

14.00000 

0.52.433 

0.70000 

16.00000 

0.59600 

0.700G0 

10.00000 

0.66333 

0.70000 

20.00000 

0.73067 

0.90000 

0.0 

0.02500 

0.90000 

2.00000 

O.099OU 

0.90000 

4.00000 

0.17300 

0.90000 

6.00000 

0.24900 

0.90000 

0.00000 

0.32500 

0.90000 

10.00000 

0.40150 

0.90000 

12.00000 

0.47800 

0.90000 

14.00000 

0.53300 

0.90000 

16.00 0 0 0 

0.62500 

0.90000 

lrt.oooun 

0.69400 

0.. 900 0 0 

2U. 00000 

0.760U0 


0.09000 

0.2 »0'ju 

n.ii 

* 

0.02106 

3.0/500 

0,0 

• * 

0.02260 

6.22750 

o . n 

• • 

0.027bl 

7.37000 

0.0 

!* « 1 

0.05486 

7.63750 

0.0 

* « » t 

0.04333 

7.50750 

0.0 

IJ. • 

0.05382 

/.2/OJO 

0.0 

n.M 

0.06713 

6.78750 

0.0 

1! , *} 

0 . 08660 

6.0 / 300 

0.0 

(K:l 

0 .lit 98 

5.19750' 

0.0 

n,5i 

0.13304 

3.95000 

0.0 

0 « J 

0.03446 

0.65000 

n .o 

It • l 

0.02341 

5.02500 

0.0 

'll . 

0.02495 

6.38250 

0,0 

0 . o 

0 . 0505? 

7.29000 

0.0 

n ♦!» 

0,05870 

7.51250 

0.0 

u . .1 

0.04892 

7.50250 

0.0 

[}«() 

0.06450 

6.57000 

0 . 0 

u.o 

0 . 08556 

5.76250 

0.0 

n . 

0.11417 

4.90500 

0.0 

U.O 

0 . 1 5296 

4.09250 

0.0 

0 m § \ 

0.21292 

3.25000 

0.0 

O.'i 

0.02649 

0.61667 

0.0 

0.0 

0.02600 

3.38333 

0.0 

o.u 

0.02891 

5.70667 

0.0 

0.0 

0.03638 

6 . 53353 

0.0 

0.0 

0.04707 

6.48333 

0.0 

u,u 

0.06254 

6.10000 

0.0 

u.o 

0.08509 

5.32000 

0.0 

0.0 

0.11566 

4.53333 

0.0 

u.o 

0.15698 

3.79667 

0.0 

0 . u 

0.27156 

2.44267 

0.0 

u.o 

0.27230 

2.68333 

0.0 

0.0 

0.03333 

0.75000 

0.0 

o.o 

0.03145 

3.15000 

0.0 

0.0 

0.03604 

4,00000 

0.0 

n.u 

0.04Q32 

5.10000 

0.0 

o.o 

0.071 43 

4.55000 

0.0 

n . o 

0.10080 

3,98000 

u.o 

u.o 

0.13580 

3.52000 

0.0 

u.o 

0.17839 

3.10U00 

0.0 

U.u 

0.22836 

2.75000 

0.0 

U - 0 

2.79037 

0.24800 

o.o 

n . '•* 

0 . 03778 

2.25000 

0.0 



Table F-l : 04 0a Aerodynamic Characteristics, Clean vehicle 

No Landing Gear, No Speed Brake 
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MACH 

ALPHA 

CL 

CD 

L/D 

S.B(%) 

L • G • 

0.30000 

0.0 

0.02250 

0.09380' 

0.23987 

0.0 

1 .001)00 

U..500UU 

2.00000 

0.08162 

0.02486 

3.28279 

u.o 

1 . U0lj:;0 

0.30000 

4.00000 

0.14075 

0.02640 

5.33116 

0.0 

1.00000 

0.30000 

b. 00000 

0*20350 

0.03141 

6.47*343 

0.0 

l.noooo 

0.30000 

8.00000 

0.26O25 

0.0386 -> 

6.88681 

0.0 

1 . Oil Mu 0 

0.30000 

lO.OCUUl) 

0.32875 

0.04713 

6.97571 

0.0 

unoimy 

o .3onno 

12.00000 

0.39125 

0.05762 

6. 791162 

0.0 

1 

0.30000 

14.00000 

0.45562 

0.07093 

6.4238 j 

0.0 

l.nor.o 

0.30000 

lb. 00000 

0*52000 

0.08940 

5.81677 

u.o 

l.orjn.jo 

0.3U00Q 

18,00000 

0,58200 

0.11578 

5.02691 

0.0 

1.09 ! : U 

0*30000 

20.00000 

0.64-iOU 

0 . 1 6884 

3.6o0()3 

0.0 

3 .09000 

0.50000 

0.0 

0.01550 

0.03824 

0.40529 

0.0 

1.0 09 90 

0.50000 

2.00000 

0.08487 

0.02721 

3.11802 

0.0 

i.ooo io 

O.bOOOO 

4.00000 

0.15425 

0.02876 

5.36532 

0.0 

1 .ooo ou 

0.50000 

6.00000 

0.22260 

0.03432 

6.48286 

0.0 

1.00900 

0.50000 

8.00000 

0.29075 

0.04250 

6.84083 

0,0 

l.OC 

0.50000 

10.00000 

0.35725 

0.05272 

6 • 7761.6 

0.0 

l.OUU-JO 

0.50000 

12.00000 

0.42375 

0.06830 

6 ■ 2 U 4 * j- () 

0.0 

1*000 HJ 

0.50000 

14.00000 

0.49107 

0.08916 

5.51690 

0.0 

1 .OOOJO 

0.50000 

lb. onooo 

0.56000 

0.11797 

4.74700 

0.0 

1 .000 JO 

U.5U0HU 

lH.OOUOU 

U . O2b00 

0 « ioo fo 

3.99330 

0.0 

i • n n c • • o 

0.50000 

20.00000 

0.69200 

0.21672 

3.19302 

0.0 

1 . 00 J'lO 

0.70000 

0.0 

0.01633 

0.03029 

0.53929 

0.0 

l.ncijoo 

0.70000 

2.00000 

0.09067 

0.03060 

2.96315 

0.0 

1.0O-*:)0 

0.70000 

4,00000 

0.16500 

0.03271 

5.04378 

0.0 

3 . oounu 

0.70000 

b. 00000 

0.23767 

0.04018 

5.91541 

0.0 

1.00090 

0.70000 

8.00000 

0.31033 

0.05167 

6.00649 

0.0 

1.0UMJO 

0.70000 

10.00000 

0.38150 

0.06634 

5.75060 

0.0 

1.00000 

0.70000 

12.00000 

0 • 452o7 

0.088 19 

5.092-57 

0.0 

l .ousno 

0.70000 

14.00000 

0.52433 

0.11946 

4.30913 

0.0 

l.OOOrO 

0.70000 

16.00000 

0.39600 

0.16078 

3.70693 

0.0 

1.010)00 

0. 7000U 

18.00000 

0.66333 

0.27536 

2.40896 

0.0 

1. onooo 

0.70000 

20.00000 

0.73067 

0.27610 

2.64640 

0.0 

1.80f' ( :0 

0.90000 

0.0 ' 

0.02500 

0.03713 

0.67325 

0.0 

1.004! 0 

0.90000 

2.00000 

0.09900 

0.03323 

2.81022 

0.0 

1.00090 

0.90000 

4.00000 

0.17300 

0 . 03984 

4,342i.9 

0.0 

1 .000 

0.90000 

6.00000 

0.24900 

0.05262 

4.73172 

0.0 

l.uo no 9 

O.9QU00 

8.00000 

0.32500 

0.07523 

4.32017 

0.0 

1.0049V 

0.90000 

10.00000 

0.40150 

0.10468 

3.83532 

0.0 

l.Ot*. 1 JU 

0.90000 

12.00000 

0.47800 

0. 139o0 

3.42418 

0.0 

1.00(100 

0.90000 

14.00000 

0.53300 

0.18219 

2.03534 

u.o 

l.ono-;o 

0 .90000 

16.00000 

0.62800 

0.23216 

2.70499 

0.0 

l.OOV'O 

0.90000 

18.00000 

0 • 6940 U 

2.80217 

0.24766 

0.0 

1 .00 ’>'-0 

0.90000 

20.00000 

0.76000 

0.34158 

2.22497 

0.0 

1 . 00 000 
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MACH 


ALPHA 


S.B (%)- 


L.G 


0.30000 

0.30000 

0.30000 

0.30000 

0.30000 

0.30000 

0.30000 

0.30000 

0.30000 

0.30000 

0.30000 

O.bOOOO 

0.30000 

0.50000 

0.50000 

0.50000 

0. 50000 

O.bOOOO 

0 . 500U 0 

O.bOOOO 

0.50000 

0.50000 

0.70000 

0.70000 

0.70000 

0.70000 

0.70000 

0.70000 

0.70000 

0.70000 

0.70000 

0.70000 

0.70000 

0.90000 

0.90000 

0.90000 

0.90000 

0.90000 

0.90000 

0.90000 

0.90000 

0.90000 

0.90000 

0.90000 


, 0.0 

2.00000 

4.00000 

. 6.00000 

8.00000 

10.00000 

12.00000 

14.00000 

16.00000 

15.00000 
20.0UU00 

0.0 

2.00000 

4.00000 

6.00000 

а. ooooo 

10.00000 

12.00000 

14.00000 

16.00000 
Id. ooooo 

20.00000 

0.0 

2.00000 

4.00000 

б . 00000 

8. 00000 

10.00000 

12.00000 

14.00000 

16.00000 

13.00000 

20.00000 

0.0 

2.00000 

4.00000 

6.00000 

8.00000 

10.00000 

12.00000 

14.00000 

16.00000 

18.00000 

20.00000 


CL 


0.00375 
0.06287 
0.12200 
0 . 1847b 
0.24750 
0.31000 
0.37250 
0 . 43b87 
0.50125 
0 . 56325 
0.62525 
- 0.00325 
0.06612 
0.13550 
0.20375 
0.27200 
0.33850 
0.40500 
0.47512 
0.54125 
0 . 6072b 
0.67325 
- 0.00242 
0.07192 
0.14625 
0.21892 
0.29158 
0.36275 
0.43392 
0.50558 
0 . 5772b 
0.64458 
0.71192 
0.00625 
0.08025 
0.15425 
0.23025 
0.30625 
0.38275 
0.45925 
0 . 53425 
0.60925 
0 . b7525 
0.74125 


CD 


0.09630 

0 .02736 

0.02890 

0.03391 

0.04116 

0.04963 

0.06012 

0.07343 

0,09190 

0 . 11828 

0.16934 

0 . 0“074 

U. 02971 

0.03123 

0.03682 

0.04500 

0.05522 

0.07080 

0.09166 

0.12047 

0.15926 

0.21922 

0.03279 

0.03310 

0.03521 

0.04268 

0.05417 

0.06884 

0.09139 

0.12196 

0.16328 

0.27786 

0.27860 

0.03963 

0.03773 

0.04234 

0.05512 

0.07773 

0.10718 

0.14210 

0.18469 

0 . 23466 

2 . 6U4r . t 

0.34408 


L/D 


0.03894 
2 . 29763 
4.22125 
5.44793 
6.01299 
6.24649 
6.19624 
5.94978 
5.45450 
4.7623 3 
3 . 69232 
- 0.07977 
2.22540 
4.33607 
5.53349 
6.04415 
6.12985 
5.72053 
5.16186 
4 . 49285 
3.83268 
3.07107 
- 0.07373 
2.17284 
4.15323 
5 . 12955 
5.38311 
5.26939 
4.74809 
4.14543 
3.53534 
2.31981 
2.55535 
0.15770 
2.12703 
3.64298 
a.17698 
3.93999 
3.57112 
3.23198 
2.89273 
2 . 59627 

0.24076 

2.15431 


25 . 000 on 

25 . 00000 ' 

25.00000 

25.00000 

2 b. 000 Lin 

25.00000 
2b« 0IJ0OO 

25.00900 

25.09000 

25.00009 

25.00000 

26.00000 
25 . 0-30 00 
25,00090 

25.00000 

25 . 00000 

25.00000 

25.00000 

26.00000 

25.00000 

25.00000 

25.00000 

25.00000 

25.00000 

25 . 00000 

25.00000 

25.00000 

25.00000 

25.00000 

25.00000 

25.00000 
25 . UOO0O 

25.00000 

25.00000 
25.00900 

25.00000 

25.00000 

25.00000 

25.00000 

25.00000 

25.00000 

25.00000 

25 .0000 (3 

25.00000 


1 1,00090 
1 1VUOOOO 
1,00000 
1 . 000.90 
1.00000 

1.00000 

] .OOOOO 

i .nnoon 

1.00000 
i. no ooo 
. i.nouoo 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 
1 . 000:30 

1.00000 

1. OOOOO 

1.00000 

1 .00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 
l.onooo 

1.00000 

1.00000 

1.00000 
l.onooo 
l.nonoo 

1. 00000 

1.00000 
l.ouonu 
l.nonoo 

1.0 oii no 

3 .00000 

1. 00000 
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MACH 


0.30900 
0.3U0UD 
0.30900 
0.30000 
0.30003 
0.30.100 
0.30100 
0.30000 
0. 30000 
0.30-iOQ 
0.30000 
0.500)0 
0.50000 
0.50000 
0.50000 
0.50000 
0.50000 
0.50000 
0.50000 
0.50000 
0.50000 
0.50000 
0.70000 
0.70000 
0.70000 
0.70000 
0.70000 
0.70000 
0.70000 
0.70000 
0.70000 
0.70000 
0.70000 
0.90000 
0.90000 
0.90000 

0.90000 
0.90000 
0.9000U 
0 .901)00 
0.90000 
o.9.j;u'io 
0.90000 
0.90000 


Table 


ALPHA 

CL 

CD 

0.0 

-0.01500 

0.09880 

2.00000 

0.04412 

0.02906 

4.00 100 

0.10325 

0.0314U 

6 . OOwO) 

0.16 '.00 

0.0564 1 

O.OUihJU 

0.22875 

0.04308 

10.00000 

0.29] 25 

0.05213 

12.00000 

0 . 55 373 

0.062*2? 

14.00000 

0.41«12 

0 . 07n93 

16.0000(1 

0.46230 

0.09440 

io.ooooo 

0.541- 50 

0.12078 

20.00 100 

0.60650 

0.17184 

0.0 

-0.02200 

0,04324 

2.00000 

0.04737 

0 . 0321*1 

4.0:1000 

0.11675 

U . 05 -73 

6.00000 

0.18500 

0.03932 

8.00090 

0.25323 

0.04750 

10.00000 

0.31973 

0.03772 

12.00000 

0.38625 

0.07550 

14.00000 

0.45437 

0.09416 

lb. 00000 

0.52250 

0.12297 

18.00000 

0.58850 

0 . It->176 

20.00000 

0.65430 

0.22172 

0.0 

-0.02117 

0.03629 

2.00000 

n n *<i7 

V « !J\K' l / 

0.03660 

4.000C0 

0. 12730 

0.03771 

6.00O0U 

0.20017 

0.04318 

8.00)00 

0.27283 

0.035 ^7 

10.00000 

0.344U0 

0.07154 

12.00300 

0.41517 

0.09369 

14.00 000 

0.48d83 

0.12446 

16.00000 

0 . 53o60 

0.16578 

18.00000 

0.62583 

0. 28036 

20.00000 

0.O9317 

0.281 10 

0.0 

-0.01250 

0.04213 

2.00300 

0.06150 

0.04023 

4.DJOOO 

0.13550 

0.04*-H4 

6.00300 

0.21150 

0.05762 

8.03000 

0.20750 

0.06023 

10.00000 

0 . 36*1 00 

0.10960 

12.00.TJ0 

0.4 a f)50 

0. 14--60 

14.00000 

0 • 51550 

0.18719 

16 , 00000 

0.59050 

0.23716 

18.00000 

0.65350 

2.80 717 

20.00030 

0. 72250 

0 . 34658 


>'L/D 

S . B (%} 

L.G. 

-0.16182 

50.000*10 

1.0--2.JO 

1.47760 

50. DO.). )0 

l.OOO'jiJ 

3.2-V.07 

50.00)00 

1 « 0-»:'*)U 

4.36.094 

bU.nO'TM) 

] .06000 

3.25924 

50 • 90*101) 

1 .(ICO 10 

5.30728 

50.00600 

1,0 Mi 00 

5.-/1942 

DO.O'iOUU 

i.o-iooo 

5.30693 

50. 00 J 00 

1 .09,(00 

5. 11 <41 

50.00-19.1 

J . oi'OOO 

4.3605J 

hO.'I'lOO 

1. , fn'Oijo 

3.32949 

50. ‘h* >■)') 

L.O.SOOO 

-0.50874 

50.000 JO 

1.00030 

1.470o4 

50. OOOOO 

1 • 09 '00 

3.43931 

50.0-U'n'j 

! .OP’J'UJ 

4.7()403 

50,«1w0o0 

1 .03 no 

5. 5. >134 

50 »*K’00n 

1 .0)0*10 

3.33932 

50. ()*»*):)') 

1.00000 

5.269-31 

50.00000 

l.OOOUO 

4 . 82807 

50.00000 

1,03000 

4.24903 

50.00000 

3 .00000 

3.63005 

50.000.j0 

1. 00 000 

2.95180 

50.00000 

J .09(300 

-0.39005 

50 .00000 

1. 00000 

1 . 49355 

50. 00300 

1.00000 

3.30075 

50 .iJ0i)00 

1.33000 

4 . 43067 

50.00 1 00 

i.noo no 

4.81 473 

50.09000 

1 .noono 

4.82192 

50.00000 

1,00000 

4.42195 

50,00000 

1. 00 000 

3.91 151 

5 0,0 O') 00 

1 .noono 

3 * 3o093 

50.00000 

1.030.30 

2.23224 

50.00000 

1.00300 

2.46592 

50.00000 

1.00000 

-0.29630 

50.09000 

1,09000 

1.52076 

50, OOOOO 

1.03930 

3.02174 

50.00000 

1 ,90 WOO 

3.67037 

50.09030 

1 .0.9000 

3.505ol 

5U . OoO 0 0 

1.00900 

3.31877 

50.0. /WOO 

l.l!) 000 

3.0*1343 

50.0-,'O30 

1. 1)9 300 

. 2.75393 

50.0. ■000 

1.0 Jl’.jO 

2.40904 

50.03600 

1.0911:10 

0.23387 

50.06000 

L. 11**900 

2 . 08467 

60.03000 

1.03*100 


F-4: 040a Aerodynamic Characteristics, Speed Brake = 50% 
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CL 


CD 


L/D 


MACH 


0-. 30000 
0.30000 
0.3UOUO 
0.30000 
0.3000'J 
0.30000 
0.30000 
0.30000 
0.300 00 
0.30000 
0.30000 
0.50000 
0.50000 
0.50000 
0.50000 
0.50000 
0.50000 
0.50000 
0.50000 
0.50000 
0.50000 
0.300 00 
0.70000 
0.70000 
0.70000 
0.70000 
0.70000 
' 0.70000 
0.7000 0 
0.70000 
0.70000 
0.70000 
0.70000 

o.yonoo 

0.90000 

0.90000 

0.90000 

0.90000 

0.90000 

0.901)00 

0.90000 

0.90000 

0.90000 

0.90000 


ALPHA 


0 .0 

2.00000 

4.00000 
fi, 00000 

8.00000 

10. no nno 

12.00000 

14.00000 

16.00000 

18.00000 

20.00000 

0.0 

2.00000 

4.00000 

6.00000 

8.00000 

10.00000 

12.00000 

14.00000 

10 .00000 
10 . 0000 ') 
20.00008 

0.0 

2.00000 

4.00000 

6.00000 

8.00000 

10.00000 

12.00000 

14.00000 

16.00000 

18.00000 

20.03000 
0.0 

2.00000 

4.00000 

6.00000 

8.00000 

10.00000 

12.03000 

14.00000 

16.00000 

18.00000 

20.00000 


— 0-. 0-3-375 
0.02337 
0.08430 
0.14725 
0.21000 
0.27280 
0 . 33500 
0.39937 
0.4637b 
0.52875 
0.58773 
-0.04075 
0 . 02362 
0.09800 
0.16425 
0.23430 
0.30103 
0.36750 
0.43562 
0.50375 
0.56975 
0.63575 
-0.03992 
0.03442 
0.10875 
0.18142 
0.25408 
0.32525 
0 • 39642 
0.46808 
0.53978 
0.60708 
0.67442 
-0.03125 
0.04275 
0.11675 
0.19275 
0.26875 
0.34525 
0.4217b 
0.49675 
0.57175 
0.637/5 
• 0.70375 


0-* 1-0-130 
0.03236 
0.03390 
0.03693 
0 . 04ol6 
0.054 63 
0.06512 
0.07043 
0.09690 
0.12388 
0.17434 
0.04574 
0.034/1 
0.03625 
0.0438? 
0.05000 
0 . 0602 ? 
0.07580 
0 . 0966o 
0.12547 
0.36426 
0.22422 
0.037/9 
0.03810 
0.04021 
0 . 04768 
0.05937 
0 . 07384 
0.09639 
0.12666 
0.16828 
0.28286 
0 . 28360 
0.044o3 
0.04273 
0.04734 
0.06012 
0.08273 
0.11218 
0.14710 
0.18969 
0.23966 
2.809o 7 
0.34908 


-0 . 333-1 7 
0.78404 
2.49252 
3.78418 
4.54930 
4 ,98' -30 
5.14*'.*8 
5.09251 
4.78603 

4. 2.o 4 7-; 
3.5/1 33 
-0.89082 
0.82450 
2.70349 
3.97520 
4.o6 u 80 
4.99821 
4.84843 
4.50682. 
4.03493 
3.46853 
2 ♦ £33*) V") 
-1.05637 
0.90317 
2.70431 
3.80507 
a.294 39 
4.40474 
4.11273 
3,68681 
3.20746 
2.14622 
2.37807 
-0.70015 
1.00050 
2.46811 
3.20590 
3.24858 
3.0// fc'J 
2 • 8o7 3 9 
2.61879 
2 . 38564 
0 . 22698 
2.01603 


S.B(%) 


70 .00 ')•).) 
75. 000')'.) 

75.000. ) J 

75.00000 
75.030*11) 
75.00500 
75 .0:80)0 
73.0-.ii , 't , 0 
75 . QOtii.i i 
7O.0'9:)J 
75.006 * .1 
75.000ir) 
75.00909 
78, 0090.1 

75.000. ' i 
7!5,OOOiH! 
75.00 '.'"’J / 
75. 0*3090 

75.00000 
75.00900 

75. OOO-in 

75,0:1000 

75.00 .-;o 
Vh.Ow r-0 

7o. OO?:. 0 

75. 00908 
75.000:10 
75. 00 non 
75.00690 
75. 00 9 DO 
75.00 600 

75.00000 
75. 00!) 09 
75.00090 

75.00000 
75.000*10 
75.00009 

75.0000. ) 

75.U9i)9i) 

75. 00 00 . 1 

75.00000 
75 . OOOrO 

75.00000 

7 5.0. HSU 0 


L.G. 


1 .'000 00 

1.0- JOOO 

l.onuoo 

1.00090 

i .oouoo 

1,00390 

1.00. 1no 
l.nnooo 

1 • Oil's jll 
1.00990 
1.09 i.K) 

1 .IJ.i.J.'O 

i.iv. i-o 

1.0 :o:jo 
1 . 03:3 oo 

1.0 OJ 90 

1.000. 8) 

1.0 1000 
1.00600 
1.00 -.00 
1.00900 

l.onooo 

1 . 00 9 »)0 
l.O’r-90 

1.00- 490 
1 • 06 JOO 
1,0».JOO 
1 . O j 390 

1.00000 
l.Lli.3 10 

l .ooooo 

1.00. 190 

1 .00. 300 

1.00 . 4 90 

1.0. j J.iO 
1,00:100 
l.OODuO 

1 .00- 600 
1 . 05990 

l.orooo 
1.06 000 
l.OOOOO 
1. CO-509 
1 « 0.1.990 


Table F-5: 040a Aerodynamic Characteristics, Speed Brake - 75% 
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MACH 


0.30000 

0. 30000 

0.30000 

0.30000 

0.30000 

0.30000 

0.30000 

0.30000 

0.30000 

0.30000 

0.30000 

0.30000 

0.30000 

0.30000 

0.30000 

0.50000 

0.30000 

0.30000 

0.30000 

0 « 30000 

o . bnoui) 

0 • 50000 

0.70000 

0.70000 

0.70000 

0.70000 

0.70000 

0.70000 

0.70000 

0.70000 

0.70000 

0.70000 

0.70000 

0.00000 

0.90000 

0.90000 

0.900'K) 

0.90000 

0.90000 

0.90000 

0.90000 

0.90: !,')() 

0 .90000 

0.900 JO 


Table 


*•’• 'o.cn if!*: 


ALPHA 

CL 

CD 

0.0 

-0.032.30 

0.10380 

2.00000 

0.0 0662 

0.0348:) 

4. 00000 

0.06575 

0.03640 

6.00000 

0.12650 

0.04141 

8.00000 

0.19123 

0.0486O 

10.00000 

0.25373 

0.05713 

12. OOO.TO 

0.51625 

0.06762 

14, 00000 

0.36062 

0.00093 

16. OOOOJ 

0.4*500 

0.099*0 

16.00000 

0.50700 

0.12378 

20.00000 

0.56900 

0 . 3 7 684 

0.0 

-0.05950 

0.0482* 

2.00000 

0.00967 

Ovj 

JO 

O 

• 

O 

4.00000 

0.07923 

0 . 03873 

6.00000 

0.14750 

0.04432 

a.ooouo 

0.21575 

0.05230 

10.0000U 

0.26225 

O.Oe.272 

12.00000 

0.3467b 

0.07830 

14.09000 

0.41687 

0.09916 

16.00001) 

0.48500 

0 . 12797 

10. onooo 

0.53100 

0. 16: -7b 

20.00000 

0.61700 

0 . 22672 

0.0 • 

—0 • 05867 

0.0*029 

2.00000 

0.01367 

0.0*060 

4.00000 

0.09000 

0.04271 

6.00000 

U. 16267 

0 . 030 16 

B.ooonu 

0.23533 

0 . 063 67 

10.00000 

0.30630 

O'. U7r>3* 

12,00000 

0.37767 

0 . 09809 

14.00000 

0.44933 

0.129*6 

16.00000 

0.52100 

0.17076 

18.00000 

0.58833 

0.2833d 

20.00000 

0.65567 

0.28630 

0.0 

-0.05000 

0.04713 

2.00000 

0.02400 

0.0*523 

4.00000 

0.09800 

0.04984 

6.00000 

0.17400 

0 . 06262 

6.00000 

0.25000 

0.00323 

10.00000 

0.32650 

0.11468 

12.00000 

0.40300 

0.14960 

14.00000 

0.47600 

0.19219 

lb. 00000 

0.55300 

0.24216 

18.0000.) 

0.61900 

2.81217 

20.00000 

0.6H300 

0.35158 


L/D 

S.B(%) 

L.G. 

-0.303/6 

) 00.000!)!* 

1.00000 

0 . l'Hu)2 

100.00000 

1 . 00000 

3 .8062.5 

too . 000 00 

1 . 00000 

3.10297 

100.00000 

1.00000 

3.9302d 

100. 000 Oil 

l.ouono 
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